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departure  and  spin  resistance  characteristics. 

The  procedure  is  based  on  nonlinear  lifting  line  theory  which  has  been  modified  to 
include  unsteady  wake  effects.  A discrete  vor'ex  lattice  representation  is  used  for  the 
time  dependent  wake,  whereas  a finite  element  representation  is  used  to  describe  the 
time  dependent  wing  load  distribution,  which  is  assumed  concentrated  along  a single 
line  at  the  25%  chord  position.  In  accordance  with  lifting  line  theory,  each  chordwise 
section  is  assumed  to  behave  like  a two-dimensional  airfoil  at  an  effective  angle  of 
attack  defined  by  geometry  and  induced  flov-  angularity.  The  usual  assumptions  of  lift- 
ing line  theory  apply  to  the  method,  viz.  moderate  to  high  wing  aspect  ratio,  moderate 
sweep  angle  and  incompressible  flow.  Additional  assumptions  are  also  implicit  regard- 
ing effects  of  the  unsteady  t’iow  or.  section  chordwise  load  distributions  and  stall. 

Arbitrary  nonlinear  section  lift  curves  may  be  introduced,  and  time  histories  of  the 
overall  wing  longitudinal  and  lateral  aerodynamic  characteristics  and  spamvise  load 
distribution  may  be  examined  as  the  wing  is  pitched  (up  or  down)  through  the  stall  range. 
Lateral  asymmetries  (e.  g.  aileron  deflection)  may  also  be  introduced  and  removed 
during  the  pitching  motion  to  simulate  effects  of  control  actuation. 

Calculations  are  presented  which  compare  the  method  with  existing  theory  and  test 
data  and  which  show  effects  of  various  planform  modifications  and  airfoil  section 
variations  on  wing  stalling  characteristics.  For  airfoils  which  have  a negative  lift 
curve  slope  (subsequent  to  stall),  the  solutions  sometimes  exhibit  abrupt  losses  in  lift, 
lift  hysteresis,  and  asymmetric  stall  (even  with  zero  lateral  asymmetry).  These 
results  are  due  to  the  existence  of  multiple  solutions  to  the  lifting  line  equations,  which 
have  been  known  to  occur  with  negative  section  lift  curves  (e.  g.  References  I - 3).  The 
transient  flow  nature  of  the  present  formulation  may  provide  a means  of  choosing 
between  the  various  possible  loadings  at  a given  angle  of  attack. 


The  theory  predicts  that  lift  hysteresis  during  stall  is  strongly  influenced  by  wing  plan- 
form  shape  and  by  the  steepness  of  the  negative  lift  curve  slope  beyond  stall,  increasing 
aspect  ratio  and  decreasing  the  section  negative  lift  curve  slope  tend  to  reduce  stall 
hysteresis.  Laterally  asymmetrical  span  load  distributions  are  produced,  even  at  zero 
yaw,  by  introducing  and  then  removing  an  aileron  deflection.  This  gives  rise  to  a ' zero 
beta"  rolling  and  yawing  moment  not  unlike  those  observed  in  tests  at  post  stall  angles 
of  attack.  The  zero  beta  asymmetries  are  shown  to  be  alleviated  by  the  same  factors 
which  reduced  stall  hysteresis.  Definitive  wind  tunnel  tests  involving  dynamic 
measurements  of  aerodynamic  force  and  moment  coefficients  and  of  span  load 
distribution  are  required  to  evaluate  the  accuracy  of  the  method. 
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ABSTRACT 


A computational  procedure  has  been  developed  for  predicting  the  time  dependent 
longitudinal  and  lateral  aerodynamic  characteristics  of  wing-body  configurations  at 
angles  of  attack  up  to  and  beyond  stall.  The  purpose  of  the  procedure  is  to  provide 
the  aircraft  designer  with  a tool  for  simulating  and  alleviating  such  adverse  wing 
stalling  characteristics  as  wing  rock,  wing  drop,  loss  of  roll  control  or  roll  control 
reversal,  etc. , and  thereby  should  lead  to  the  design  of  aircraft  with  improved  stall, 
departure  and  spin  resistance  characteristics. 

The  procedure  is  based  on  nod  inear  lifting  line  theory  which  has  been  modified  to 
include  unsteady  wake  effects.  A discrete  vortex  lattice  representation  is  used  for 
the  time  dependent  wake,  whereas  a finite  element  representation  is  used  to  describe 
the  time  dependent  wing  load  distribution,  which  is  assumed  concentrated  along  a 
single  line  at  the  25%  chord  position.  In  accordance  with  lifting  line  theory,  each 
chord  wise  section  is  assumed  to  behave  like  a two-dimensional  airfoil  at  an 
effective  angle  of  attack  defined  by  geometry  and  induced  flow  angularity.  The  usual 
assumptions  of  lifting  line  theory  apply  to  the  method,  viz.  moderate  to  high  wing 
aspect  ratio  , moderate  sweep  angle  and  incompressible  flow.  Additional 
assumptions  are  also  implicit  regarding  effects  of  the  unsteady  flow  on  section  chord- 
wise  load  distributions  i*nd  stall. 

Arbitrary  nonlinear  section  lift  curves  may  be  Introduced,  ard  time  histories  of  the 
overall  wing  longitudinal  and  lateral  aerodynamic  characteristics  and  spanwlse  load 
distribution  may  be  examined  as  the  wing  is  pitched  (up  or  down)  through  the  stall 
range.  Lateral  asymmetries  (e.  g.  aileron  deflection)  may  also  be  introduced  and 
removed  during  the  pitching  motion  to  simulate  effects  of  control  actuation. 

Calculations  are  presented  which  compare  the  method  with  existing  theory  and  test 
data  and  which  show  effects  of  various  planform  modifications  and  airfoil  section 
variations  on  wing  stalling  characteristics.  For  airfoils  which  have  a negative  lift 
curve  slope  (subsequent  to  stall),  the  solutions  sometimes  exhibit  abrupt  losses  in 
lift,  lift  hysteresis  , and  asymmetric  stall  (even  with  zero  lateral  asymmetry'). 

These  results  are  due  to  the  existence  of  multiple  solutions  to  the  lifting  line  equations, 
which  lave  been  known  to  occur  with  negative  section  lift  curves  (e.  g.  References 
1-3).  The  transient  flow  nature  of  the  present  formulation  may  provide  a means  of 
choosing  between  the  various  possible  loadings  at  a given  angle  of  *t.'iack. 

The  theory  predicts  that  lift  hysteresis  taring  stall  Is  strongly  lifluenced  by  wing 
planform  shape  and  by  the  steepness  of  the  negative  lift  curve  slope  beyond  stall. 
Increasing  aspect  ratio  and  decreasing  the  section  negative  lift  curve  slope  tend  to 
reduce  stall  hysteresis.  Laterally  asymmetrical  span  load  distributions  are  produced. 
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even  at  zero  yaw,  by  Introducing  and  then  removing  an  aileron  deflection.  This 
gives  rise  to  a "zero  beta"  rolling  and  yawing  moment  iut  unlike  thoue  observed 
In  tests  at  post  stall  angles  of  attack.  The  zero  beta  asymmetries  are  shown  to  be 
alleviated  by  the  same  factors  which  reduced  stall  hysteresis.  Definitive  wind  tunnel 
tests  Involving  dynamic  measurements  of  aerodynamic  force  and  moment  coefficients 
and  of  span  load  distribution  are  required  to  evaluate  the  accuracy  of  the  method. 
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INTRODUCTION 

As  Is  well  known,  there  are  a large  number  of  aerodynamic  effects  which  may  limit 
the  maneuverability  of  a fighter,  especially  as  thrust  to  weight  ratios,  structural 
placards  and  physiological  limits  on  the  pilot  are  raised.  These  include  such  phenom- 
ena as  shock  induced  buffet,  wing  rock  and  wing  drop  (probably  due  to  asymmetrical 
separation  with  a subsequent  rolling-yawing  motion),  nose  slicing  (possibly  due  to  the 
establishment  of  asymmetrical  flow  and  yawing  moments  in  the  nose  region),  pitch-up 
Goss  of  longitudinal  stability  due  to  flow  separation  at  the  tips),  loss  of  lateral  stabil- 
ity at  high  angles  of  attack,  loss  of  roll  control  and  roll  control  reversal  at  high  angles 
of  attack,  and  large  drag  increases  (which  lower  the  maximum  sustained  load  factor). 

It  is  generally  believed  that  some  type  of  flow  separation  is  responsible  for  each  of 
these  phenomena.  The  present  report  deals  with  the  development  of  a simplified  aero- 
dynamic tool  for  estimating  the  aerodynamic  characteristics  of  wings  and  wing-body 
configurations  of  moderate  to  high  aspect  ratio  undergoing  transient  pitching  motion 
with  flow  separation.  Hence,  the  method  developed  in  this  study  should  prove  useful 
for  the  prediction,  simulation  and  possible  inhibition  of  those  adverse  stalling 
phenomena  included  within  the  framework  of  the  theory. 

As  noted  above,  the  present  formulation  Is  for  wings  of  moderate  to  high  aspect  ratio 
to  which  lifting  line  theory  may  be  applied.  In  this  regard  there  has  appeared  some 
interesting  analytical  work  (1-3),  based  on  Prandtl  lifting  line  theory,  which  indicates 
that  multiple  solutions  of  the  lifting  line  equations  are  possible  under  certain  conditions. 
The  multiple  solutions  are  predicted  only  if  the  two-dimensional  lift  curves  have  either 
discontinuities  or  regions  of  high  negative  slope,  beyond  stall,  as  illustrated  in  Figure 
1.  One  interesting  feature  of  the  multiple  solutions  is  that  they  can  yield  span  load 
distributions  which  are  either  symmetrical  or  asymmetrical  about  the  wing  centerline 
even  at  zero  sidelsip  angle.  According  to  Sears  ^),  the  asymmetrical  solutions 
suggested  to  von  Karman  that  large  rolling  moments  could  be  produced  near  the  stall 
without  postulating  any  initial  rolling  velocity  (or  other  asymmetries  such  as  slight 
sideslip  or  wing  panel  misalignment).  Sears  states  that  the  usual  "textbook"  explana- 
tion for  the  large  rolling  moments  near  stall,  which  is  based  on  an  initial  rolling 
velocity  which  stalls  one  panel  and  destalls  the  other,  fails  to  account  for  the  violence 
of  the  rolling  moments  experienced  in  a wind  tunnel  with  a model  held  fixed  to  the  sting. 
The  establishment  of  any  one  particular  load  distribution  (either  sy nun etr leal  cr 
asymmetrical)  probably  depends  on  the  relative  stability  of  the  corresponding  circulation 
distribution  to  small  disturbances. 

The  present  report  attempts  to  apply  nonlinear  lifting  line  theory  to  the  modeling  and 
prediction  of  aerodynamic  effects  during  stall  penetration  at  high  angles  of  attack. 
Previous  discussions  of  stall  aerodynamics  have  noted  the  occurrence  of  asymmetrical 
flow  separation  at  zero  yaw  angle  as  a possible  forcing  mechanism  for  wing  rock  and 
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and  other  types  of  pre  and  post-stall  motions  (e.g.  References  4-7).  It  Is  postdated 
herein  that  the  time  dependent  zero  beta  yawing  and  rolling  moments  predicted  by  the 
nonlinear  lifting  line  theory  are  contributory  to  the  limit  cycle  type  of  motions  apt  to 
occur  during  stall.  This,  of  course,  will  require  careful  experimental  verification. 

The  lifting  line  formulation  utilized  herein  is  a finite  element,  unsteady  wake,  incom- 
pressible flow  theory  and  is  somewhat  more  general  than  the  General  Dynamics/Convair 
nonlinear  lifting  line  procedure  reported  in  References  8 - 10.  Besides  excluding  time 
dependent  wake  effects,  the  later  procedure  was  limited  to  laterally  symmetric  load 
distributions  at  zero  yaw  angle,  did  not  include  effects  of  yaw,  and  could  i»t  handle 
body  interference  effects.  Nevertheless,  the  procedure  was  found  to  correlate 
reasonably  well  with  wind  tunnel  data  through  the  stall  (Ref.  9). 

The  procedure  as  formulated  herein  may  be  applied  at  either  zero  or  non  zero  yaw 
angle  (yaw  is  handled  by  skewing  the  wing  panels).  Both  symmetrical  and  asymmetrical 
load  distributions  may  be  obtained,  even  at  zero  yaw  angle.  Three  dimensional  unsteady 
aerodynamic  effects  are  included  by  allowing  shed  vortices  in  the  wake  to  vary  in 
strength  with  distance  and  time.  The  strengths  of  the  shed  vortices  are  related  to 
those  of  the  corresponding  bound  elements  at  an  earlier  time,  based  on  the  convective 
time  dalay  at  free  stream  velocity'  between  the  bound  vortex  and  the  particular  wake 
station.  Although  the  theory  is  unsteady  from  the  point  of  view  of  wake  induced  effects, 
it  is  assumed  that  the  two-dimensional  airfoil  chordwise  loadings  and  sectional 
characteristics  in  stall  are  steady  state.  Thus,  the  assumption  is  implicit  that  such 
two  dimensional  phenomena  as  associated  with  the  dynamic  stall  of  helicopter  blades 
are  of  a shorter  time  scale  than  the  unsteady  wake  effects  (see  also  Section  2.  2). 

The  present  approach  makes  use  of  the  key  assumption  in  lifting  line  theory,  viz. 
that  each  choi'dwise  section  acts  like  a two-dimensional  airfoil  at  an  effective  angle 
of  attack  equal  to  the  local  geometric  angle  of  attack  less  the  induced  angle  of  attack. 

In  order  for  this  assumption  to  remain  valid,  (L,  e.  for  the  chordwise  load  distribution 
to  remain  two-dimensional)  the  aspect  ratio  must  be  sufficiently  large  to  limit  span- 
wise  flow  effects  , and  vortex  effects  associated  with  high  leading  edge  sweep 
angles  must  not  be  dominant  in  the  stall.  Besides  the  assumptions  noted  above  the 
method  requires  the  selection  of  (or  switching  between)  the  various  multiple  span  load 
solutions  which  may  exist  at  a given  angle  of  attack.  Some  of  these  solutions  may  be 
indicative  of  unstalled  load  distributions,  some  of  partly  stalled  distributions  and  some 
of  nearly  fully  stalled  distributions.  The  criterion  used  for  selecting  among  the  various 
matematically  possible  solutions  is  empirical,  and  it  is  recognized  that  this  aspect  of 
the  formulation  requires  further  analysis.  Because  of  this  empiricism,  and  as  a 
result  of  the  assumptions  noted  previously,  the  need  for  wind  tunnel  tests  involving 
dynamic  measurements  of  aerodynamic  forces,  moments  and  span  load  distribution 
can  not  be  over-emphasized. 

The  mathematical  basis  for  the  present  theory  is  formulated  in  Section  II.  Results  of 
some  simplified  analyses  and  properties  of  the  solutions  are  also  presented  in  Section  II. 
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The  method  of  solution  and  a dieussion  of  the  stability  and  convergence  of  the  iteration 
procedure  are  presented  in  Section  in.  More  extensive  numerical  results  involving 
comparison  with  existing  test  data  and  linear  theory,  are  found  in  Section  IV.  Also 
presented  in  Section  IV  are  parametric  studies  showing  effects  of  planform  shape 
and  of  two-dimensional  lift  curve  characteristics  on  stalling  properties.  Mathemati- 
cal details  of  the  formulation  and  a listing  of  the  computer  program,  along  with 
detailed  instructions  for  its  use,  may  be  found  in  Reference  12. 


* The  reader  is  also  referred  to  Reference  25,  which  came  to  the  attention  of  the 
author  just  prior  to  completion  of  the  present  report,  and  which  also  applies  non- 
linear lifting  line  theory  to  the  design  of  wings  which  minimize  aircraft  departure 
and  spin  entry  tendencies. 
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THEORY  AND  ASSUMPTIONS 


The  present  section  presents  the  mathematical  basis  and  assumptions  for  the  nonlinear 
lifting  line  theory  with  a time  dependent  wake  as  utilized  herein.  The  theory  may  be 
used  to  obtain  the  time  dependent  span  load  distribution  and  corresponding  time 
dependent  wing  force  and  moment  coefficients  for  wings  and  wing-bodies  of  arbitrary 
planform  experiencing  a predetermined  schedule  of  pitching  motion  and/or  of  lateral 
control  deflections. 

The  present  lifting  line  theory  is  similar  to  that  used  by  Weissinger  in  his  L-method^1! 
but  has  been  modified  to  include  nonlinear  section  aerodynamics,  body  interference 
effects  and  time  dependent  wake  effects.  Effects  of  sideslip,  or  yaw,  are  represented 
as  differential  right  and  left  sweep  angles  (like  a skewed  wing).  No  fuselage  cross 
flow  terms  are  included  in  this  sideslip  representation.  The  theory’  as  presented 
below  allows  for  the  arbitrary  chordwise  positioning  of  the  wing  control  points,  although 
the  bound  vortex  location  is  always  retained  at  the  25%  chord  location. 

2. 1 VORTEX  SYSTEM 


The  vortex  system  used  in  the  analysis  is  pictured  in  Figures  2 and  3,  and  is  assumed 
planar  except  for  the  vertical  displacement  of  the  image  vortices  in  the  fuselage.  The 
exposed  wing,  as  shown  in  Figure  2,  is  segmented  into  N equal-span  elements  (up  to 
20)  of  width  h in  the  spanwise  y direction.  Each  element  is  a parallelogram  in  shape 
and  consists  of  a bound  vortex  segment  along  the  1/4  chord  line  (which  may  be  swept), 
two  streamwise  trailing  segments  and  a closing  shed  vortex  element  at  a downstream 
distance  ix-  Vo  it,  where  V®  is  the  free  stream  velocity  and  At  is  the  time  step. 
In  the  time  dependent  solution  the  strength  F-  ^(t)  of  this  vortex  element  represents 
the  circulation  carried  by  the  wing  at  the  j 'tll’spanwise  element  (j  - 1, 2 ...  N)  at  the 
current  time  t.  The  downstream  trailing  vortex  elements  are  of  identical  shape, 
but  are  of  varying  strength  as 


(t) 


r 

‘j,  k-i 


( t - At) 


(1) 


where  k is  an  index  describing  the  vortex  element  number  in  the  wake  (see  Figure  3). 

If  the  maximum  number  of  downstream  parallelogram  elements  to  be  included  is  M-] , 
then 

1 £ k £ M 

where  k = 1 corresponds  to  the  wing  element  at  the  current  time  t,  and  k = M refers 

Preceding  page  blank 
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to  the  final  wake  element,  which  is  assumed  unclosed.  If  M = 1,  the  solution  reduces 
to  that  for  a steady  state  wake. 

Representation  of  the  fuselage  (assumed  an  infinite  circular  cylinder)  by  Image  vortex 
elements,  as  depicted  In  Figures  3 and  4,  has  found  extensive  use  in  the  literature, 
e.g.  References  13-16.  The  body  Images  are  located  at  identical  axial  positions  as 
the  primary  (wing  and  wake)  vortices.  Since  only  the  trailing  vortex  segments  are 
properly  imaged  (cancellation  of  primary  vortex  induced  velocities  normal  to  the 
fuselage)  the  lateral  and  vertical  locations  of  the  image  vortices  is  based  on  location 
of  the  primary  trailing  vortices.  As  shown  in  Figure  4,  the  images  are  located  along 
a straight  line  in  the  y-z  plane  between  the  primary  trailing  vortex  and  the  fuselage 
center  at  a radial  distance  Rj  from  the  centerline 

KI  * kb2/rp  (2> 


where  Rg  is  the  body  radius  and  R~  is  the  radius  to  the  primary  vortex. 

The  velocities  induced  by  the  connecting  bound  elements  between  the  image  vortices 
are  also  included,  based  on  Ref.  13,  even  though  the  bound  image  vortices  do  not 
satisfy  the  no  flow  condition  through  the  fuselage.  However,  these  velocities  were 
not  turther  augmented  by  the  factor  suggested  in  Reference  14,  since  it  was  believed 
that  the  justification  for  this  correction  was  insufficient-  The  bound  image  vortices 
were  included  in  the  evaluation  of  the  overall  force  on  the  wing-bodv  combination,  in 
accordance  with  the  Trefftz-plane  relationships  given  in  Reference  16. 

2.2  LIFTING  LINE  EQUATIONS 


The  lifting  line  assumption  states  that  each  spanwise  station  i of  a wing  acts  like  a 
two-dimensional  airfoil  at  an  effective  angle  of  attack  equal  to  the  local 

geometric  pitch  angle  of  attack  (includes  effects  of  wing  incidence,  twist,  control 
deflection,  and  body  cross  flow)  less  the  downwash  angle  of  attack  ag.  induced  by 
the  trailing  vortex  systems,  viz. 


a 

EFF^ 


(3) 


Thus,  for  an  airfoil  with  lift,  pitching  moment,  nnd  profile  drag  curves  of  functional 
form 


CL  * CL  <V'  CM  " °M  CD  - ° 


D <"„>  ’ 
P 


(4) 


lifting  line  theory  assumes  that 


6 


(5) 


CL  (t) 

= °L. 

a 

EFF. 

(t) 

i 

1 

i 

J 

CM  » 

II 

o 

S 

^EFF, 

(t)l 

1 

i 

L t 

J 

cD  (t)  = Cr 
P. 

i 


D r^EFF  (t)] 

psL  i J 


Equations  (5)  Indicate  also  that  the  unsteady  effects  to  be  Included  in  the  theory  enter 
only  through  the  time  dependence  of  <*£FFt  * viz*  through  changes  with  time  in 
c*p.  and/or  a£y  As  noted  in  the  Introduction,  the  basic  section  characteristics  as 
denned  functionally  by  Equations  {4)  are  assumed  independent  of  time.  This  requires 
that  various  linear  and  nonlinear  two-dimensional  unsteady  flow  effects,  such  as 
dynamic  stall  (e.g.  Reference  17),  be  of  an  appreciably  shorter  time  scale  than  the 
wake  and  geometric  effects  Included  herein.  Thi3  assumption  is  equivalent  to  assuming 
that  the  chordwise  pressure  distribution,  which  affects  the  nature  of  flow  separation 
on  the  airfoil,  rapidly  approach  the  steady  state  distribution  for  the  current  angle 
of  attack.  Combining  unsteady  two-dimensional  effects  with  the  nonlinear  lifting 
line  theory  was  assumed  beyond  the  scope  of  the  present  effort. 

The  control  point  for  evaluating  oeff.  Is  taken  at  an  axial  distance  xcp.  from  the 
leading  edge.  In  che  usual  Prandtl  lifting  line  formulation  x^.  - c./4  .l  Herce, 
there  would  be  no  induced  velocity  contribution  fx*om  the  wing  bound  vortex  segments. 

In  order  to  generalize  the  formulation  for  arbitrary  x,^,  the  downwash  angle 
<*£)  has  been  defined  as  the  difference  between  the  total  three-dimensional  downwash 
angle  ^30^  at  the  control  point  and  an  equivalent  two-dimensional  downwash  angle 
<*20-  from  an  infinite  span  bound  vortex  along  the  c/4  line  and  equal  in  strength  to 
Fj  lit).  Thus 


o > = a -a 
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In  terms  of  downwash  velocity  components  Aw  induced  by  the  parallelogram  lattice 
elements  of  unit  circulation  strength 
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where  for  each  j,  k lattice  element  and  at  each  control  point  i 
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The  subscripts  in  Equation  (8)  are  T trailing  element,  B bound  or  shed  element,  I 
Image,  L left,  R right,  F forward  and  A aft.  Mathematical  expressions  for  each  of 
the  elements  in  Equation  (8),  based  on  the  law  of  Biot  and  Savart,  are  given  in 
Reference  12. 


The  equivalent  two-dimensional  dowmva^h  angle  is 


to2Dtri,l  <« 


(t) 


* VVA  C0B(Af  W 


(9) 


where  is  the  quarter  chord  sweep  angle  (assumed  negative  for  the  left  wing 
panel.  1 £ i £ N/2,  and  positive  for  the  right  wing  panel,  N/2  + 1 s.  I s N),  and  8 
is  the  side  slip  angle.  The  usual  sign  convention  taking  the  downwash  as  positive 
has  been  employed  in  the  above  expressions. 


The  pitch  angle  of  attack  is  composed  of  the  elements 

“p  (t)  = <*(t)+<*T  + AorR  (t)  + AorB  (t) 
i i i i 


(10) 


where  a (t)  is  the  angle  of  attack  of  the  root  section,  »T.  is  ike  built  in 
geometric  twist,  \a  (t)  is  the  effective  roll  asymmetry  angle  of  attack  due  to 
deflection  of  ailerons,  and  Ac^  (t)  is  the  additional  angle  of  attack  due  to  the 
body  upwash.  i 
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Here  iw  is  the  wing  root  incidence  angle,  and  the  remaining  quantities  are  defined 
in  Figure  4. 

The  unknowns  are  the  current  time  values  of  tne  circulation  on  the  wing  jit),  and 
may  be  found  by  solving  Equation  (3)  and  the  first  of  Equations  (5).  Thus,  inserting 
Equal  *'ons  (6),  (7),  (9)  and  (10)  into  Equation  (3)  gives 
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Expressing  Cj^.  (t)  in  Equation  (5)  in  terms  of  the  bound  circulation  gives 


2 Fi  l(t) 

V c,  = CL  ^EFFj  W 3 (13) 
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Equation  (13)  implies  that  the  Kutta-  Joukowskl  law  Lj  (t)  - pVu>.r  j (t)  holds  in  the 
unsteady  as  well  as  in  the  steadv  state  case.  This  assumption  is  discussed  further 
in  Section  2.5. 

Equations  (12)  and  (13)  constitute  the  set  of  lifting  line  equations  to  be  solved  for 
^EFF-  W M i^)-  The  k (t)»  1,  are  known  from  solutions  at  previous 

times  through  Equation  (1).  iiefore  describing  the  general  method  of  solution 
several  simplified  solutions  of  these  equations,  which  do  not  require  any  computer 
calculations,  are  discussed  below. 

2.3  SIMPLIFICATIONS  FOR  A LINEAR  LIFT  CURVE 

In  general.  Equations  (12)  and  (13)  must  be  solved  by  an  iterative  method,  because 
of  the  nonlinear  nature  of  the  lift  curve  as  expressed  in  Equation  (13).  It  is 
precisely  these  nonlinearities  which  lead  to  the  occurancc  of  multiple  span  loadings, 
lift  hysteresis,  and  zero  beta  rolling  and  yawing  moments.  However,  it  is  beneficial 
to  examine  several  simplifications  which  occur  when  the  lift  curve  is  linear,  in 
order  to  compare  the  theory  with  existing  methods  and  to  better  understand  the 
nature  of  the  complications  due  to  the  nonlinear  effects. 
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2.3.1  LINEAR  LIFT  CURVE  - In  the  case  of  a linear  lift  curve  Equation  (13) 
becomes 


(t) 
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where  a^  is  the  section  lift  curve  slope.  Combining  Equation  (14)  with  Equation  (13) 
then  gives  a single  set  of  linear  equations  which  may  be  solved  directly  for  F.  (t). 

L«  1 


The  result  is 
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where  the  repeated  indices  indicate  a summation,  and  6..  is  the  Kronecker  delta 
symbol . lJ 


Equation  (15)  represents  a linear  set  of  algebraic  equations  which  may  be  solved 
directly  for  each  step  in  time.  If  the  problem  is  initiated  from  a steady  state 
condition,  then  M = 1 for  t = 0 and  the  double  summation  term  vanishes.  In  matrix 
notation  Equation  (15)  now  becomes 


[AW){r]  (0)  J -{ap  (0)) 


06) 


where  elements  of  the  [AW]  square  matrix  are  given  by  the  term  in  square  brackets 
in  Equation  (15),  whereas  j Fj  (0)}  and  {ap^  (0)}  are  column  matrices  made  up  of 
the  elements  r.  j (0)  and  ap.  (0)  respectively.  Inverting  Equation  (16)  gives 

fri  (0)}=  (AWj'hofp  (0)}  (17) 

At  subsequent  time  steps  terms  in  the  double  summation  plays  a role  in  the  solution. 
Thus,  according  to  Equation  (1),  with  t = At,  k = 2 now  corresponds  to  the  solution 
T.  , (0).  Hence,  we  require  M a 2 in  order  to  account  for  unsteady  wake  effects. 

2.3.2  EQUIVALENCE  WITH  WEISSINGER  THEORY 

2.3.2. 1 Steady  State  - The  pvjaeirt  lifting  line  formulation  becomes  identical  to  the 
Weissinger  L-method^^  for  a sinear  lift  curve  in  the  steady  state  case  provided 
the  control  point  xcp  is  placed  ai  distance  equal  to  50%  of  chord  from  the  lifting  line 
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(75%  chord  for  zero  sweep  angle).  This  may  be  shown  by  considering  a single  horseshoe 
vortex  element  (k  = 1 only),  with  a bound  segment  and  two  trailing  segments.  By 
Equation  (6) 


of  = a - a 
D.  3D.  2D 

l i i 


and  from  Equation  (9)  with  (xcp  - xc/4).  cos  (A  - 0)  = (l/2)  c^ 

r 


O' 


2D.  irVa3cl 


According  to  Equation  (3) 
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o — cv  . - cx  + 

EFF.  P.  3D.  7i  V_c 


i 


(18) 


For  thin  airfoil  aerodynamics,  a.  = 2 77;  therefore  by  the  steady  state  form  of  Equation 
(14) 


_LL 


= O' 


IT  V c.  EFF. 

CO  l I 


and  substitution  in  Equation  (18)  yields 


(19) 


which  is  the  flow  tangency  boundary  condition  at  the  75%  chord  control  point  used  in  the 
Weissinger  approach.  The  nonlinear  lifting  line  formulation,  with  the  control  point 
location  such  that 


(x  - x ,.)  cos  (A  - 0)  = c/2 
cp  c/4 

may  therefore  be  taken  equivalent  to  a nonlinear  Weissinger  formulation. 


(20) 


2. 3. 2.2  Unsteady  Theory  (Two-Dimensional).  It  should  be  noted  that  the  theory 
also  reduces  to  an  unsteady  form  of  the  Weissinger  theory,  both  for  two-dimensional 
and  three-dimensional  flows.  Thus,  for  a two-dimensional  unsteady  flow  Equations  (6) 
and  (7)  become 

01  — -,(t)  =~~  aw,,  r,_  (t)  - 


r1  (t) 


EFF' 


k k 


2*V»  (xcp"xc/4) 
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where  Aw^  is  the  downwaah  at  x^  due  only  to  forward  and  aft  spanwise  segments  of 
the  k'th  element.  Thus,  from  Equation  (8) 


Aw,  = Aw  + Aw 

k 13  * 


o-pjt)  = Of  (t)  - — Aw.  r.  (t)  + 1 

EFF  p k k 2^V  (X  -x 

°°  cp  c/4 

Using  the  linearized  form  of  the  unsteady  Kutta-Joukowski  law.  Equation  (14), 
assuming  the  lift  curve  slope  a = 2n,  and  taking  xC))-  x c j = c/2  then  gives 


— Aw  r (t)  = cv  (t)  (22) 

V°°  k k p k = 1, 2, . . . 

Equation  22  states  that  the  instantaneous  downwash  angle  at  the  75 r,o  chord  position 
from  all  bound  and  shed  spanwise  vortex  segments  is  equal  to  the  instantaneous 
angle  of  attack. 

Equation  (22)  may  be  expressed  in  terms  of  the  differential  vortex  strengths. 
Referring  to  Figure  5 


Thus,  from  Equation  (21) 


AVk<t>“(4w: b + 4wb  )rk(t) 

Fk  \ 

and  since  Aw  = -Aw 

JD  . 

A,  F, 

k k+l 


substitution  in  Equation  (22)  yields 


- >»y>rk(t,=Vt, 


whe/°  as  seen  from  Figure  5, 
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ai'j  (t)  - rx  (t)  - r0(t)  = rj(t) 


is  the  instantaneous  bound  vortex  strength  on  the  airfoil  at  the  25%  chord  position. 
As  also  noted  in  Figure  5,  the  force  free  convection  of  the  wake  requires 


Af  k (t  + At>  Ark-1  (t) 


(24) 


whereas  the  condition  that  the  total  vorticity  (bound  and  shed)  remain  zero  for  all 
time  requires 

Af  (t  + At)  = - [ryt  + At)  - fj(t)  ] (25) 


Equation  (25)  requires  that  the  first  time  step  be  taken  with  k = 2,  for  which  the 
first  shed  vortex  is  at  a distance  e = ( n - 3/4)  c from  the  trailing  edge,  as  shown 
in  Figure  5.  Here  n is  the  number  of  chord  lengths  per  time  step. 

Equations  (22)  and  (23)  both  show  that  the  unsteady  theory  reduces  to  satisfying  a flow 
tangency  condition  in  the  two-dimensioral  case  when  the  control  station  is  taken  at  the 
75%  chord  position. 

The  equivalence  of  the  present  unsteady  wake  lifting  line  formulation  with  the  classical 
solution  of  Wagner  (Reference  18)  for  a two  dimension  airfoil  in  unsteady  flow  is 
readily  shown  through  use  of  Equations  (23)  - (25).  Thus,  Figure  6 compares  the 
calculated  value  of  the  change  in  lift  coefficient  AC^(t)  due  to  a step  change  in 
angle  of  attack  A^p,  as  calculated  from  Equation  (23)  - (25),  with  the  Wagner  function. 
Lift  coefficient  was  calculated  from  the  bound  circulation  strength  fj  (t)  through  the 
unsteady  Kutta-J oukowski  law.  Equation  (13). 


The  step  size  used  in  the  compution  for  the  shedding  of  discrete  vortices  was  one 
chord  length.  The  ciose  correspondence  between  the  discrete  vortex  formulation 
and  the  continuous  vortex  sheet  theory  of  Wagner  shows  that  very  little  accuracy 
is  lost  in  the  discrete  vortex  model.  It  ?s  also  clear  that  unsteady  wake  effects,  due 
to  the  variation  of  shed  as  well  as  of  trailing  vorticity  are  accounted  for  in  the  formu- 
lation. Effects  of  wake  vorticity  on  choidwise  load  distribution  are  not  included  as 
was  mentioned  previously.  Incorporation  of  these  effects  would  require  addition 
bound  vortex  elements  and  control  point  locations  on  the  wing. 

2. 3.  2. 3 Unsteady  Theory  (Three-Dimensional) 

The  unsteady  finite  element  representation  may  also  be  readily  shown  to  lead  to  a 
flow  tangency  condition  at  the  75%  chord  position  in  the  three  dimensional  case.  The 


13 


The  demonstration  is  similar  to  that  in  Section  2.  3.2. 1 for  the  steady  state  case, 
except  that  the  unsteady  form  of  the  Kutta-  Joukowski  law  is  used  as  given  by  Equation 
(13)  and  the  summation  is  carried  out  over  k = 1,  2, . . . M rows  of  vortex  elements 
in  the  wake.  The  result  is, 

“3D  (t)  = ~V~  4Wi  Ak(t)=Qp.(t)  <26) 

i 1:0  j,k  ‘i 

provided  the  control  point  is  located  according  to  Equation  (20)  and  provided  the 
steady  state  thin  airfoil  theory  lift  curve  slope  value  is  used  . 

2.4  SIMPLIFIED  SOLUTIONS  FOR  A NONLINEAR  LIFT  CURVE 

Several  calculations  have  been  performed  with  N = 2 (one  vortex  per  wing  panel) 
and  a tri-linear  2D  lift  curve.  The  two-vortex  representation  is  the  simplest 
possible  model  for  obtaining  asymmetrical  load  distributions  at  zero  yaw  angle 
and  should  be  illustrative  of  the  type  of  results  obtainable  with  a larger  number  of 
horseshoe  vortex  elements.  Only  simplified  steady  state  solutions  are  discussed 
below,  since  the  corresponding  simplified  unsteady  solutions  did  not  include  effects 
of  the  shed  vortex  segments.  Unsteady  nonlinear  solutions  will  be  presented  in 
connection  with  the  iterative  procedure  described  for  solving  Equations  (12)  and  (13) 
in  Section  3. 

As  illustrated  geometrically  in  Figure  7,  several  different  steady  state  load  distributions, 
both  symmetrical  and  asymmetrical,  may  be  possible  for  a wing  with  atri-iinear  lift 
curve  at  geometric  angle  of  attack  o?p  when  ofp>  ff9>  Thus,  solutions  1 through 
6 represent  the  following  loadings  (see  also  Figure  lb): 

1.  Symmetrical,  both  panels  unstalled  with  |x>sitive  lift  curve  slope  a. 

2.  Symmetrical,  both  panels  stalled  with  negative  lift  curve  slope  - b. 

3.  Symmetrical,  both  panels  stalled  with  zero  lift  curve  slope 

4.  Asymmetrical,  one  panel  unstalled  and  other  panel  stalled  with  lift  curve 
slope  -b. 

5.  Asymmetrical,  one  panel  unstalled  and  other  panel  stalled  with  lift  curve 
slope  0. 

6.  Asymmetrical,  both  panels  stalled,  one  with  lift  curve  slope  -b  and  the 
other  with  lift  curve  slope  0. 
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The  determination  or  which  load  distribution  will  occur  depends  on  the  airfoil  2-D 
lift  characteristics  (b/a,  ofj  and  a 2),  on  the  geometric  angle  of  attack,  and  on  wing 
aspect  ratio. 


Symmetrical  (unstalled  with  lift  curve  slope  a) 

For  loading  1,  the  lifting  line  equations,  [Equations  (12)  and  (13)  ],  simplify, 
respectively,  to 


C = a {«  - y C -y  C ) 
L ' P rll  L /12  L } 

1 1 Cl 


(27) 


provided  the  effective  angle  of  attack  a , which  is 

EFr 


a = a - y C -y  C 
EFF  p yil  L y12  hn 
v 1 2 


(28) 


is  less  than  o^.  Here  y^  is  the  downwash  angle  induced  at  control  point  1 (assumed 
at  the  center  of  bound  element  1)  due  to  horseshoe  vortex  1 of  strength  such  that 
= 1.  Similarly,  ?l2  is  the  downwash  angle  induced  at  control  point  1 due  to  horse- 
shoe vortex  2 of  strength  such  that  Ct  = 1.  The  influence  coefficients  ^ are  relate 
to  the  downwash  terms  Aw^/V^  used  m previous  sections  through  the  expression 
Tij  = (Va>Cj_/2)  Awy/V*.  Here  also,  the  control  point  position  has  been  taken  on  the 
bound  vortex,  so  that  Aw2p  in  Equation  (12)  cancels  the  contribution  from  the  bound 
vortex. 

For  symmetrical  loadings,  Ci,  = C ^ . From  the  Biot-Savart  law  the  coefficients 
yjj  and  y12  are  readily  shown  to  be 


c /c 
yil  ” vJR 


_ V! 

yi2  3tt  JR 


where  and  c2  are  the  respective  chords  and  c is  the  average  chord.  • For  a wing 
of  symmetrical  planform  with  only  two  horseshoe  vortices  c^/c  = c2/c  = 1. 

Solving  Equation  (27)  for  /ace^  gives 


for  the  unstalled  3-D  lift  curve.  Equation  (29)  may  be  compared  with  the  analogous 
lifting  line  result  with  a constant  downwash  angle  for  a - 2 tt.  The  latter  Is 


ap  , 

C^J 2tt  = - + 2/j^  whereas  Equation  (29)  gives 


C /2tt  

L 1 4-  — — — 

3 m 


^mmetrical  (stalled  with  lift  curve  slope 


For  loading  2,  the  effective  angle  of  attack  defined  by  Equation  (28)  must  lie  between 
anda^*  The  lifting  line  equation  is  now 

°L1  = a‘Vb  VT11  CL  -ri2CL ,-“l> 

1 ' i.  Ct 

Solving  for  C /a  a gives,  assuming  symmetry, 

L1  1 


1 + ~ (1  -or  /a A 
a p i 


aa.  1 - (b/a)  ay^  (1  + y^/y^) 


Symmetrical  (stalled  with  lift  curve  slope  0) 

For  loading  3,  the  three  dimensional  lift  coefficient  is  simply  given  by 

CT/aai  = 1 - ACT/aa, 

LI  L i 


where  AC^is  the  drop  in  lift  coefficient  between  ando^.  For  loading  3 we  require 
thatoEpp.  as  given  by  Equation  (28)  exceed 

Asymmetrical  (panel  #1  unstalled,  panel  #2  stalled  with  dCT  /do>  = -b) 

For  loading  4,  the  lifting  line  equations  give 


C = a (a  - y.  „ C - y „ C \ 

Lx  *p  'll  1^  r12  L > 


C,  = aa,  - b (a  - y„  C - y C - a ) 

L 1 *p  21  L,  y22  L0  r 

^ 12 
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where 


a - a - y C - y C 
EFFX  p rll  Lx  T12  L2 


(33) 


a = a - 7 C - 7 C 
EFF2  P 21  L1  *22  L2 

We  require  that  *ess  al  that  aEFF2  ^ie  between  a^  anda2.  Solving 

Equations  (32)  and  (33)  with 

y = 7 and  7 = 7 due  to  the  symmetrical  wing  planform  yields 

±U  ul  11  UCd 


c /aa  = N /D, 
L1  1 1 


(34) 


C /aa  * N /D 
L2 


and  C /C  = N./N 
L2  L1  1 


(35) 

(36) 


where 


/«p\  r b 

n = Ur)  i -- 
1 r 1/  I & 


(aTn-aV 


-d*-)  ar12 


I,avn-aV-(1  + I> 


and  D = ( 1 + a rlx)  (1  - J a Tjj)  + ~ (a  V^)2 


+ (!  + -)  (1  ♦at'jj) 


Equation  (36)  has  been  plotted  in  Figure  8 for  = 4.  Also  indicated  in  Figure  8 is 
the  maximum  angle  of  attack  limit  above  which  aEEE  exceeds  ap  As  seen  from 
Figure  8,  Equation  (36)  predicts  an  asymmetrical  rolling  moment  at  angles  of  attack 
below  that  for  3-D  when  b/a  s 1.5.  Hence,  some  deterioration  in  aircraft  de- 

parture characteristics  might  be  expected  under  these  conditions. 
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Asymmetrical  (panel  #1  unstalled,  panel  #2  stalled  with  dC^/do-p  = 0) 
For  loading  #5,  the  lifting  line  equations  reduce  to 


(37) 


We  also  require  that  ttgpp.  be  less  than  a 1 and  that  appp2  be  greater  than«2*  where 
aEFFi  a°d  «£pp9  are  defined  according  to  Equations  (33).  Solving  Equations  (37) 
gives 


AC. 


\ c V**!1  ~ (Via/Yll)  * T11  ('  ~ I^T 


acr 


1 + a y 
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AC. 


a«. 


= 1 - 


acr 


CL  /CL 
L2  L1 


AC. 


-1 


(Y12/ril>aYll 


-1 


_1 

1 + ay 
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(38) 


(39) 


(40) 


Equation  (40)  has  been  plotted  in  Figure  9 for  values  of  AG^/actj  between  0 and  1 and 
/R  = 4.  Also  shown  in  Figure  9 is  the  maximum  angle  of  attack  at  which  G'gpjr 
<a  (unstalled  panel  limit)  and  the  minimum  angle  of  attack  at  which  a ppp0  > o ^ 
(stalled  panel  limit).  The  latter  limit  is  dependent  upon  b/'a.  Figure  9 sho\s  that 
with  b/'a  > 1.5,  asymmetric  loadings  of  type  5 may  exist  over  an  a zone  which 
increases  in  width  with  increasing  b/a  and  AC  /a a, . P 

.Li  j. 

Asymmetrical  (panel  #1  stalled  with  dC^/da  -b,  panel  #2  stalled  with  dCj/da  = 0) 


In  this  case  the  lifting  line  equations  give 


C = a«t  - b (cp  - T12  C - yn  C - 
1 {*  \ 


(41) 


C - aa  - ACL 
L2  1 
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with  the  requirement  that  oppp,  lie  between  and  g2»  and  that  appp2  be  greater 
thanof^.  Equations  (41)  were  solved  for  C^/C^  giving 


% 1 

iI+a 

r ru  V AC_ 

I 1 . + ay  (l  - — -) 

[ ai  \i  11  av. 

V 

l(- 

Ac%\ 
aai  / 

» «* 
1-(I)aYu 

In  summary,  we  have  demonstrated  the  possibility  of  six  different  steady  state  load- 
ings near  stall  for  wings  with  a tri-linear  section  lift  curve.  Taking  a wing  with 
AR  = 4,  b/a  - 4 and  AC^/aaj  = 0. 6,  it  is  readily  shown  that  all  six  loadings  will  satis- 
fy the  lifting  line  equations  at  an  angle  of  attack Qp  - 1,3  o^,  which  is  just  below  the 
value  «p  = 4/3  for  the  maximum  3-D  lift  coefficient.  The  values  of  Cj^  and  C l2 
for  these  six  loadings  are  listed  in  the  table  on  Figure  10.  The  corresponding  values 
of  o;pppi  and«ppp0  are  shown  In  the  plot. 

The  question  naturally  arises;  which  of  the  six  possible  loadings  will  occur  physically  ? 
This  question  was  posed  by  Sears W who  stated  "The  choice  between  the  various  sol- 
utions must  involve  the  question  of  the  relative  stabilities  of  the  flows. " The  question 
will  be  readdressed  in  connection  with  the  iterative  method  of  solving  the  time  depen- 
dent nonlinear  lifting  line  equations. 

2.5  GENERALIZED  FORM  OF  THE  KUTTA-JOUKOWSKI  LAW  FOR  UNSTEADY 
FLOW 

A generalized  form  of  the  two-dimensional  Kutta-Joukowski  law,  viz. 

Lfy  = Polity  (43) 

has  been  used  throughout  the  formulation  in  Sections  2.2  and  2. 3 and  is  the  basis  for 
Equation  (13).  In  the  present  section  we  derive  Equation  (43)  from  the  condition  that 
the  negative  rate  of  change  of  total  momentum  associated  with  the  spanwise  vortaK 
segments  gives  the  magnitude  of  the  lift.  The  total  momentum  perpendicular  to  x 
for  a two-dimensional  vortex  system  composed  of  discrete  elements  AI^  (t)  spaced 
along  the  x axis  at  distances  x^  from  the  leading  edgec  as  shown  in  Figure  5b,  is 

PArk(t)xk 

where  the  double  index  implies  a summation.  Hence,  the  unsteady  two-dimensional 
lift  becomes 
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nt)  = -p  ~ [Ark (t) Xk  ] 


(44) 


where  the  subscript  l has  been  dropped  for  simplicity  of  rotation  [ also  Eq.  (1) 
of  Reference  (1)  ]. 

The  conditions  that  the  shed  vortex  elements,  k s 2,  are  convected  downstream  at 

M-l 

velocity  Va  , and  that  the  total  vorticity  JjT  Ark  (t)  remain  zero  for  all  time  are  given  oy 

k^l 

Equations  (24)  and  (25),  respectively.  The  instantaneous  bound  vortex  strength  is 
Aryt)  = r1  (t). 

Referring  to  Figure  5b,  taking  time  steps  with  n = 1 (one  chord  length  between  ele- 
ments), placing  the  bound  vortex  at  the  25%  chord  position,  the  momentum  at  time 
t = cAffl  is 

pAiyc/VJ  = Ar  (c/VJ  c/4  + Af  (c/4)  5 c/4 
= - r1  (cA’j  c 

where  we  have  made  use  of  Equation  (25). 

Similarly,  for  t = 2c /V^ 

pArk(2c/VJ  3^  = Ari(2cAeo)c/4  + AF2(2 !c/VJ  5c/4  + Af^c/VJ  9 c/4 

= -V1(2c/Vjc-ri  (e/VJ  c 
where  we  have  made  use  of  Equations  (24)  and  (25), 

For  t = (M  - 1)  c/V  , the  result  generalizes  to 

CO 

M-] 

PAI^[(M-1)  c/V^]  x^  = -c  ^ Tj  (kc/VJ  (45) 

k=l 

Taking  t - {3/l-2)c/V  and  subtracting  the  result  from  Equation  (45)  gives 

p | Afk  [M-Dc/Vj  - Afk  [(M-2)c/Vw]  xk  } = -cr^M-Dc/Vj 
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Dividing  by  At  = c/V^,  we  obtain 


pA{Arkt(M“i)  c/v»]xk} 

At 


-vra  ri  [(m-i)c/v  J 


and  Equation  (44)  in  difference  form  becomes 


L [(M-l)cA  ] = pV  r,[(M-l)c/V  ] (46) 

CO  CD  1 00 

which  is  identical  to  Equation  (43).  The  same  result  may  be  obtained  for  arbitrary 
step  sizes  n,  and  in  the  limit  n - 0 the  generalized  Kutta- Joukowski  law  is  tound. 
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METHOD  OF  SOLUTION  AND  ITERATION  PROCEDURE 


The  set  of  nonlinear  lifting  line  expressions  as  given  by  Equations  (12)  and  (13) 
constitute  a system  of  2N  algebraic  equations  in  the  2N  unknowns  T{(t)  and  "gpp.  (t), 

1 = 1,2,.,.,  N,  which  must  be  solved  for  each  step  in  time.  The  equations  are  nonlinear; 
because  the  functional  relation  between  ryt)  and  cvg it),  expressing  the  two- 
dimensional  lift  curve  at  station  i,  is  generally  nonlinear  in  the  stall  region  and 
beyond.  Hence  an  iteration  procedure  is  required  for  their  solution. 


The  nonlinear  lifting  line  equations  and  corresponding  wing  and  fuselage  geometric 
relations  were  programmed  for  solution  on  the  CDC  Cyber  70  co mputsr.  Details  of 
the  computer  program  and  of  the  method  of  solution  have  been  presented  in  Reference 
12,  as  was  noted  previously.  However,  the  method  of  solution  will  also  be  reviewed 
below,  because  of  the  limited  distribution  of  Reference  12  and  because  several  of  the 
assumptions  in  the  iteration  procedure  may  require  further  discussion. 

3. 1 ITERATION  PROCEDURE 

The  iteration  procedure  assumes  the  aircraft  wing  and  fuselage  geometry,  sideslip 
angle,  etc. , are  given,  and  that  the  schedule  of  pitch  angle  Op(t)  and  lateral  control 
deflection  £ceR(t)  have  been  specified.  Two  Iteration  loops  will  be  described,  one  for 
t = 0 and  the  second  for  t>  0. 

3. 1. 1 INITIAL  SOLUTION.  The  method  is  started  by  determining  the  steady  state 
solution  for  t = 0.  The  root  solving  procedure  consists  of  a simple  iteration  loop  on 
the  induced  angle  of  attack  cm  (0).  A guess  is  first  made  for  the  spanwise  variation 

of  induced  angle  of  attack  <*d^(0).  This,  together  with  the  known  (0), 
establishes  ^eff-^O)  through  Equation  (3).  The  wing  bound  vortex  strengths 
fl  (0)  are  then  found  by  table  look-up  from  the  input  aerodynamic  lift  curves 
expressed  by  Equation  (13).  For  t = 0 we  may  take  M = 1,  or  equivalently  take 
k (0)  = Tj  j (0),  hence  no  wake  vortices  need  be  considered.  The  induced  i_y’es 
are  then  recalculated,  based  on  the  downwash  velocity  components  Aw^.  k and 
A W2d.»  aJ-1  on  the  initial  iterate  for  the  vortex  strengths  ^^(O),  and  compared 
with  the  assumed  values  of  (0).  The  differences  may  be  designated  Acq)^  i°)* 

The  values  to  be  used  In  the  next  Iteration  are 

«d(2)(0)  =.or  (1)(0)  + CAor  (1)(0)  (47) 

l i l 

where  C is  a weighting  factor.  The  value  of  C affects  the  stability  and  convergence 
of  the  iteration  procedure.  Increasing  C speeds  up  convergence  (reduces  the  number 
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Precedes  page  blank 


of  required  Iterations),  but  may  destablize  the  iteration  procedure.  In  general, 
the  maximum  value  for  C depends  on  wing  aspect  ratio,  the  number  N of  spanwise 
elements,  and  on  lift  curve  slope.  Some  further  discussion  of  effects  of  C on  stability 
are  contained  in  Section  3. 1.3. 

With  an  ^ (0)  given  by  Equation  (47)  the  procedure  is  repeated,  viz. 
ut 

(2)  (2)  (2)  (2) 

aD  (0)  ^EFF  (0)  *ri  1 (0)  -*AaD  (0)’ 

i i ’ i 

until  either  convergence  is  obtained,  or  until  a maximum  number  of  iterations  has 
occurred.  Convergence  is  assumed  when  Act^  (0)  so.l  degree  for  all  i. 

L 

3. 1.2  SOLUTION  AT  SUBSEQUENT  TIME  STEPS  - The  iteration  procedure  at 
succeeding  time  steps  t>  0 is  similar  to  that  at  t = 0 with  the  following  exceptions: 


(i)  The  vortex  strengths  k (t)  the  wake  are  no  longer  identical  to  those 
on  the  wing  (k  = 1),  but  are  found  by  indexing  from  the  previous  time 
step  according  to  Equation  (1). 


(ii)  A special  logic  is  used  to  specify  the  initial  guess 


a 


(1) 


D: 


(t). 


The  later  requirement  was  found  necessary,  especially  near  stall,  because  of  the 

possibility  of  multiple  solutions,  and  because  the  iteration  procedure  apparently 

tends  to  disallow  roots  on  the  steep  negative  slope  regions  of  the  section  lift  curves. 

Therefore-,  the  wing  elements  appear  in  the  solutions  as  either  unstalled  or  fully 

stalled  (region  where  the  post  stall  lift  curve  is  near  zero).  A dominant  form  for  the 

solution  in  the  stall  region  appears  to  be  a spanwise  alternating-  pattern  of  unstalled 

ano  telly  stalled  elements.  The  significance  of  this  saw-tooth  type  of  solution  was 

not  entirely  clear.  For  example,  if  a very  large  number  of  spanwise  elements  were 

modeled,  then  such  a pattern  over  a given  spanwise  portion  of  the  wing  may  wail 

correctly  represent  partial  stall  or  heavy  buffeting  in  that  region.  However,  since 

the  present  computer  program  is  limited  to  10  elements  per  wing  panel  (Nmax  - 20), 

this  type  of  saw-tooth  stall  pattern  doss  not  appear  to  be  physically  realistic.  It 

was  found  that  the  saw-tooth  stall  pattern  could  be  avoided  in  most  cases  by  choosing 

the  initial  guess  for  a n (t)  according  to  the  following  special  logic. 

u. 

a)  When  no  wing  sections  are  stalled  at  the  prior  time  step  (i,  e.  all  converged 
values  of  ®efF[  ^ ~ At)  are  less  than  corresponding  to  the  section 

angle  of  attack  for  maximum  C^),  then  the  cq).^  (t)  are  taken  equal  to 
the  converged  values  from  the  previous  time  step  cvp.  (t  - At).  This  will 
produce  solutions  with  all  spanwise  elements  unstalled,  if  such  solutions 
exist. 
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b)  Wh  n tlie  iteration  procedure  at  the  end  rf  the  previous  time  step,  or  during 
the  Current  Iteration,  gives  a solution  with  one  or  more  stalled  sections  on 
a wing  panel  (e.g.  <*efFj  ft  " At)  eireeeds  Q^ax),  then  ti  e iteration  procedure 
is  started  (or  restarted)  with  the  assumption  that  the  induced  angles  are  zero 
on  that  particular  wing  panel.  This  procedure  tends  to  force  solutions  with 
as  many  stalled  elements  on  that  panel  as  the  lifting  line  equations  will 
allow. 

The  special  logic  represented  by  a)  forces  the  spanwise  wing  elements  to  remain 
unstalled  until  such  time  that  one  or  more  elements  must  be  stalled  to  satisfy  the 
governing  equations.  Partially  stalled  solutions  with  either  saw-tooth  stall  patterns 
or  uniform  sta’l  patterns,  are  avoided  during  this  time,  even  though  they  satisfy  the 
lifting  line  equations.  The  special  logic  given  by  b)  assumes  that  once  any  spanwise 
element  has  stalled,  it  tends  to  induce  all  adjacent  elements  of  the  same  panel  to 
stall  within  the  limits  of  the  lifting  line  equations.  Shielding  of  stall  progression 
from  one  side  of  the  airplane  to  the  other  is  assumed  provided  by  the  presence  of 
the  fuselage  (at  least  for  low  and  mid- wing  configurations).  Unstalled  solutions, 
or  solutions  with  saw-tooth  stall  patterns,  are  avoided  during  this  time  even  though 
they  are  mathematically  acceptable. 

It  should  be  noted  that  the  above  logic  will  tend  to  maximize  the  extent  of  any  stall 
hysteresis  loop  which  occurs  for  a wing  undergoing  a pitching  motion  through  the  stall. 
Thus,  stall  will  be  delayed  to  the  maximum  possible  angle  of  attack  while  o'p(t) 
is  increasing,  and  destall  will  be  delayed  to  the  minimum  possible  angle  of  attack 
while  tfp  (t)  is  decreasing. 

Thus  the  sample  computations  presented  in  Section  4 may  well  tend  to  exaggerate 
the  stall  hysteresis  loops  (and  the  angle  of  attack  range  for  the  oecurrance  of  zero 
beta  rolling  and  yawing  moments).  The  user  is  cautioned  that  the  physical  validity 
of  these  assumptions  has  yet  to  be  established,  since  there  are  at  present  insufficient 
wind  tunnel  test  data  with  dynamic  type  measurements.  An  alternative  type  of 
mathematical  formulation,  wherein  the  selection  between  the  various  multiple 
solutions  is  based  on  their  physical  stability,  rather  than  on  a special  logic  for 
the  initial  iterate,  is  clearly  desirable. 

3. 1. 3 STABILITY  OF  ITEBATION  PROCEDURE 

The  stability  of  the  iteration  procedure  may  be  readily  shown  to  be  dependent  upon 
the  spanwise  exteni  of  the  wing  vortex  elements,  upon  lift  curve  slope  a and  upon 
the  weighting  factor  C.  We  consider  for  example  the  highly  idealized  case  of  a 
wing  with  a single  vortex  element  and  no  unsteady  wake  effects  (M  - N = 1). 

Neglecting  the  subscripts  i,  j,k  the  first  iterate  at  time  t is 
/i\  Aw  * 

aD  ft)  = 7“  r ft  " At) 

V 00 
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where 


Aw 

V* 


V c 

CO 

~ 


Vco 


V 

CO 


i and 


r (t  - At)  = 


2r(t~  ad 

V e 


Is  the  normalized  circulation  (lift  coefficient)  from  the  previous  time  step.  According 
to  Sections  3. 1. 1 and  3. 1.2,  with  the  section  unstalled,  the  initial  iterate  for  r*^(t) 
is 


A (1) 

r At)  =a 


(A w \ ■'  1 


(48) 


where  a is  the  unstalled  lift  curve  slope,  and  Clq  Is  the  section  lift  coefficient  at  a =0. 
The  value  for  Aod^  (t)  which  appears  in  Equation  (47)  is  simply 


Introducing  the  weighing  factor  C into  Equation  (47)  gives  for  the  second  iteration 
in  downwash  angle 


Of 


-d(2)  (t)  - (i-C)  r(t-  at)  + c 


*r(1)  (t) 


* (2) 

The  second  iteration  for  f (t)  follows  as 


M2)  I /■>  / Aw\ ^ (1)  . 

r (t)  - a |a„(t)  - (i-c)r  (t  - At)  - c RT  r u (t)J  + cLr 


which  may  be  written,  upon  making  use  of  Equation  (48),  as 

r (2)  (t)  = (1-  o r(1)(t)  + c ja  [ ap(t)  - (yj  r (1)  (t)  ] + cLq 


Similarly,  for  the  third  Iteration 


A (3) 


(t)  = (i-c)  r (2)  (t)  + c ja  [ap(t)  - ('vr)r (2)  (t)  j + cLo 


(49) 


(50) 
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Subtracting  Equations  (49)  and  (50)  gives 


a.  (3) 

L _ tt)„ 


* (2) 
r * ; (t) 


A (2)  A (1) 

r ' ;(t)  - r w (t) 


1 - C - a C 


(51) 


The  iteration  procedure  should  be  stable  when  the  absolute  value  of  this  ratio  is 
below  one,  but  convergence  becomes  very  slow  as  the  ratio  approaches  zero. 

Since  a>  0,  the  ratio  is  generally  negative,  and  convergence  and  stability  require 


1 


Convergence  Limit 


< C< 


£ 3 


Stability  Limit 


(52) 


The  product  of  h/c  and  2v  (^w/VJ  has  been  plotted  in  Figure  3 3 versus  h/c 
for  a horseshoe  vortex  of  span  h and  with  control  point  locations  at  either  c/4  or 
3 c/4.  Examination  of  Figure  11  shows  that  the  allowable  value  of  C must  Increase 
as  the  number  of  spanwise  elements  N decreases,  or  as  aspect  ratio  increases. 

The  value  of  C should  also  increase  with  decreasing  positive  lift  curve  slope  a,  and 
C should  be  taken  somewhat  smaller  when  the  control  point  is  at  3c/4  instead  of  at 
c/4. 

Equation  (51)  was  derived  in  order  to  Illustrate  the  functional  dependence  of  C on 
paneling  geometry  and  on  lift  curve  slope,  and  is  quantitatively  correct  only  for  the 
highly  simplified  formulation  M = N = 1.  Several  trial  values  for  C mav  be  required 
in  order  to  find  the  value  which  gives  the  best  convergence  properties  for  a particular 
nonlinear  calculation.  According  to  Reference  12,  C should  range  from  about  0.1 
for  values  of  AR  =-  5 up  to  C =-0.3  for  AR  =- 12.  These  values  were  based  on  N = 15. 

A default  value  C = 0. 1 is  used  in  the  computer  program.' 

3.  2 COMPUTATION  OF  FORCE  AND  MOMENT  COEFFICIENTS 

Once  the  converged  vortex  strength  distribution  has  been  determined  at  a particular 
time  step,  total  force  and  moment  coefficients  for  the  wing  body  combination  may 
be  found  by  integrating  the  station  force  and  moment  coefficients  over  the  N wing 
elements.  The  integration  includes  contributions  from  the  bound  image  elements 
in  the  fuselage  in  order  to  account  for  fuselage  lift,  as  noted  previously.  These 
fuselage  terms  are,  however,  omitted  from  the  computation  of  drag,  side  force, 
and  rolling  moment  coefficient,  since  the  fuselage  is  believed  to  make  only  a 
negligibly  small  contribution  to  these  coefficients. 
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The  overall  coefficients  of  lift  C^,  drag  Cq,  side  force  Cy,  pitching  moment 
Cm.  rolling  moment  C and  yawing  moment  Cj^  are,  in  terms  of  the  section 
coefficients  (subscript  i), 


COS  (*£). 


M 


(53) 

(54) 

(55) 

(56) 

(57) 

(58) 


The  coefficient  ai'e  referenced  to  wing  area  S,  mean  aerodynamic  chord  c and  wing 
spaD  b.  The  angle  y refers  to  wing  dihedral.  The  effective  do wnwash  angle  ax>. 
is  taken  positive  in  the  downward  direction,  and  therefore  is  of  opposite  sense 
to  the  angle  j in  Reference  12.  Several  higher  ordered  terms,  such  as 
sin  have  been  neglected  in  the  lift  and  rolling  moment  expressions.  Furthermore, 
there  is  no  geometric  contribution  of  dihedral  angle  nor  of  sweep  angle  to  effective 
angle  of  attack.  Such  effects  are  importa  it  at  high  angles  of  attack  and  sideslip,  and 
corrections  may  be  required  to  correlate  with  test  data,  especially  with  sideslip. 
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4 

SAMPLE  CALCULATIONS 


An  extensive  seiles  of  computations  were  made  with  the  nonlinear  lifting  line  computer 
program  described  in  Sections  2 and  3.  The  computations  were  made  to  compare  with 
test  data,  with  previous  theories  and  to  assess  effects  of  wing  planform  geometry  and 
section  lift  curve  shape  on  stalling  characteristics. 

4. 1 COMPARISON  WITH  AVAILABLE  THEORY 

Several  computations  were  made  in  ordei  to  compare  the  current  method  with  existing 
theories  in  order  to  evaluate  the  accuracy  and  generality  of  the  procedure.  A 
comparison  of  the  accuracy  of  the  finite  element  representation  with  the  essentially 
exact  potential  theory  results  of  Wagner  {Reference  18)  has  been  given  previously 
in  Figure  6 for  the  case  of  a two-dimensional  airfoil  with  a linear  lift  curve  subject 
to  a step  change  in  angle  of  attack.  Because  of  the  simplifications  introduced  by 
two-dimensional  flow,  results  by  the  present  method  in  Figure  6 were  hand  computed 
instead  of  being  run  on  the  digital  computer.  The  similarity  of  the  present  finite 
element  method  with  the  exact  results  is  apparent  in  Figure  6.  Both  methods  give 
an  initial  lift  increment  equal  to  l/2  of  the  steady  state  value. 

The  present  method  is  compared  with  calculations  by  W.  P.  Jones  (Reference  19) 
and  Djojodihardjo  and  Widnall  (Reference  20)  in  Figure  12,  for  the  response  to  a 
step  change  in  angle  of  attack  for  a rectangular  wing  of  aspect  ratio  6.  The  lift  curve 
was  again  linear,  and  as  can  be  seen  from  the  figure  good  correlation  with  the  more 
exact  theories  was  found.  Both  the  present  and  more  exact  methods  show  that  the 
initial  lift  increment  is  greater  than  l/2  of  ihe  steady  state  value  when  the  aspect  ratio 
is  finite.  The  computations  in  Figure  12  were  carried  out  on  the  digital  computer  with 
values  M = 20  and  N = 14. 

Figure  13  depicts  the  variation  in  span  load  distrib-  .ion  with  time  from  the  current 
method  for  the  same  sample  computation.  Also  shown  in  Figure  13  are  the  section 
lift  coefficients  at  the  12  l/2%  and  87  l/2%  semi  span  stations  as  obtained  by  the 
method  of  Reference  20.  Once  again,  good  correlation  is  shown  with  the  more  exact 
computations. 

The  present  method  was  also  used  to  determine  the  development  of  lift  following  a 
step  change  in  angle  of  attack  for  rectangular  wings  of  various  aspect  ratios,  agaL' 
with  linear  lift  curves.  The  results  are  plotted  in  Figure  14  and  compare  well  with 
similar  variations  obtained  by  Jones  (Ref.  19)  for  wings  of  elliptic  planform.  Once 
again  we  note  that  the  starting  circulatory  lift,  increment  exceeds  one-half  the 
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steady  state  value  for  the  finite  aspect  ratio  cases,  with  the  excess  increasing  for  the 
lower  aspect  ratios.  The  computations  in  Figure  14  are  with  M = 20  and  N = 14. 

Normalized  steady  state  lift  increments  are  plotted  versus  aspect  ratio  in  Figure  15, 
to  examine  the  aspect  ratio  range  over  which  the  current  formulation  may  be  expected 
to  hold.  The  current  method  appears  to  provide  an  accurate  representation  for  lift 
down  to  aspect  ratios  below  2,  for  an  unswept  wing  of  rectangular  planform,  when 
compared  with  References  21  and  22  and  with  slender  body  theory.  However,  agree- 
ment to  such  low  aspect  ratios  is  much  less  likely  for  planforms  of  moderate  to  high 
sweep  angle  and  for  more  sensitive  coefficients  such  as  induced  drag,  section  Cl, 
etc.. 


4.2  COMPARISON  WITH  TEST  DATA 

The  wing  planform  of  the  T-2C  Navy  trainer  aircraft  served  as  the  basic  geometry 
for  the  remaining  computations  (with  various  modifications  and  adjustments).  Lift 
and  rolling  moment  data  from  Reference  23  will  be  compared  with  the  T-2C  compu- 
tations using  a steady  state  wake  (M=l),  since  the  testing  was  done  with  the  model 
held  fixed  at  a given  angle  of  attack  (pitch  and  pause). 

The  basic  T-2C  configuration  has  a wing  of  aspect  ratio  AS  = 5.07,  taper  ratio  t = 

0.495,  a quarter  chord  sweep  angle  A = 2. 27  degrees,  a wing  root  incidence  tugle 
i = 1. 7 degrees  with  2.5  degrees  of  washout  at  the  tips,  and  a dihedral  angle  T = 3°. 

The  airfoil  is  a NACA  64iA212  section  with  an  a - 0.8  (mod)  camber  line  (flaps  and 
ailerons  rigged  3 degrees  up),  and  is  constant  across  the  span.  The  wing  was 
mounted  at  mid-height  to  a fuselage  whose  shape  was  approximated  in  the  computations 
by  a circular  cylinder  of  radius  Rg=0. 119  (b/2).  The  wing  paneling  used  for  the 
computations  is  shown  in  Figure  16  in  the  case  of  N - 14.  The  number  of  spanwise 
elements  was  varied  for  different  computations.  Nonlinear  aerodynamic  section  data 
used  for  the  T-2C  computations  (and  for  the  parametric  computations)  are  plotted 
in  Figures  17  and  18,  respectively.  Curves  labeled  #1,  #2  and  #3  in  Figure  17  are 
for  a NACA  64iA212  section  at  Reynolds  number  3x10®,  6 x 1 0®  and  > 6 x 106. 

Curves  labeled  #11,  #12  and  #13  in  Figure  17  were  adjusted  for  the  -3°  trailing  edge 
rigging  of  the  T-2C  model,  and  are  to  be  used  in  comparing  with  the  wind  tunnel  data. 
Curves  labeled  #4  and  #5  in  Figure  18  are  for  a NACA  64-209  section  at  Reynolds  numbers 
3 x 10®  and  6 x 10®,  respectively,  and  were  used  in  the  parametric  calculations  together 
with  curves  #6- #8  to  examine  effects  of  the  steepness  of  the  negative  lift  curve  slope 
-b  in  the  post  stall  region.  The  airfoil  section  characteristics  in  Figures  17  and  18  were 
extrapolated  from  published  NACA  data  (e.g.  Reference  24). 

A comparison  between  the  calculated  and  measured  lift  values  for  the  T-2C  at  zero 
and  ten  degrees  yaw  angle  is  shown  in  Figure  19.  The  tests  were  conducted  tail  on, 
and  corrections  were  made  for  the  estimated  tail  lift  contribution  as  shown  in  Figure 
19b.  The  test  data  are  static  measurements  and  average  out  any  fluctuating  values 
of  Cl-  The  model  was  not  pitched  in  a manner  to  produce  lift  hysteresis  (by  taking 
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transient  measurements  with  Increasing  and  then  decreasing  with  time),  whereas 
the  calculations  were  made  with  a pitch  rate  of  8 degrees/second  {first  Increasing  and 
then  decreasing).  Even  though  It  is  difficult  to  evaluate  the  accuracy  of  the  unsteady 
method  on  the  basis  of  the  comparison  with  the  static  test  data,  the  calculations  appear 
to  indicate  a hysteresis  loop  which  extends  over  a somewhat  broader  angle  of  attack  band 
(ll°<ap<  18°)  than  might.be  expected.  This  would  be  in  accordance  with  the  special 
computational  logic  employed  in  the  method,  which  as  discussed  in  Section  3. 1. 2 may 
well  tend  to  lead  to  the  maximum  allowable  hysteresis. 

On  the  other  hand,  Reference  23  shows  the  sudden  occurence  of  a rolling  moment  at 
zero  yaw  angle  at  angles  of  attack  beyond  12  degrees.  This  may  be  taken  as  an  indi- 
cation of  the  stall  (or  partial  stall)  of  one  wing  panel,  and  may  correspond  to  the  lower 
limit  of  any  zone  of  lift  hysteresis.  The  T-2C  zero  beta  rolling  moment  data  has 
been  plotted  in  Figure  20.  Also  shown  in  Figure  20  are  computations  made  with  the 
current  method,  wherein  a roll  asymmetry  is  introduced  for  a very  short  period 
and  then  removed,  in  order  to  perturb  the  asymmetrical  loadings.  As  seen  from 
Figure  20,  the  resulting  rolling  moment  coefficients  (with  the  roll  asymmetry 
removed)  correlate  well  with  the  test  data. 

The  present  nonlinear  lifting  line  method  predicts  the  occurance  of  lift  hysteresis 
and  zero  beta  rolling  and  yawing  moments  for  a symmetrical  model  with  wing 
section  characteristics  similar  to  the  T-2C  aircraft.  Its  proper  evaluation  requires 
dynamic  wind  tunnel  testing  wherein  the  model  may  be  pitched  (up  or  down)  at  a given 
rate,  and  for  which  instantaneous  (as  well  as  time  averaged)  data  are  recorded.  It  is 
recommended  that  any  such  test  be  made  with  a wing  whose  two-dimensional  section 
characteristics  are  known.  In  addition,  a program  of  this  type  should  include 
measurements  of  span  load  distribution,  in  order  to  afford  a more  explicit  evaluation 
o£  the  theory. 

4.3  PARAMETIPC  STUDIES 

A series  of  parametric  calculations  have  been  made,  using  the  current  method,  to 
determine  effects  of  wing  planform  (aspect  ratio,  taper  ratio,  sweep  back  angle,  and 
twist)  and  of  airfoil  section  characteristics  (steepness  of  negative  lift  slope  beyond 
stall)  on  the  occurence  of  lift  hysteresis  and  zero-beta  rolling  moment.  The  purpose 
of  the  calculations  was  to  obtain  trends  and  guidelines  which  might  lead  to  the  design 
of  safer  aircraft  which  are  less  subject  to  adverse  departure  and  spin  entry  charac- 
teristics during  high  performance  maneuvers. 

4,3. 1 EFFECT  OF  NEGATIVE  LIFT  CURVE  SLOPE  AND  ASPECT  RATIO. 

Effects  of  the  abruptness  of  the  negative  lift  curve  slope  beyond  stall  and  of  aspect 
ratio  for  zero  yaw  angle  are  shown  In  Figures  21-23.  Each  plot  represents  a pitch 
hysteresis  computation  through  the  stall  with  the  angle  of  attack  first  increasing  and 
then  decreasing  at  a rate  of  8 degrees/second.  Points  with  increasing  angle  of  attack 
are  derated  by  circles,  whereas  triangles  are  used  for  points  with  decreasing  angle 


of  attack.  Where  no  triangles  appear,  the  points  are  coincident  with  the  increasing 
Op  points.  This  notation  has  been  used  throughout  the  present  section. 

Figure  21  is  for  A1  = 5. 07,  and  shows  that  the  zone  over  which  lift  hysteresis  occurs 
gradually  decreases  as  the  steepness  of  the  negative  lift  curve  slope  diminishes.  The 
same  trend  is  shown  in  Figure  22  (JR  = 8)  and  Figure  23  (At  = 12).  The  figures  also 
show  that  increasing  aspect  ratio  reduces  the  zone  of  hysteresis  for  a given  negative 
lift  curve  slope. 

A similar  set  of  pitch  sweeps  were  made  for  a yaw  angle  & = 10°.  Lift  hysteresis 
characteristics  similar  to  those  in  Figures  21-23  were  obtained.  The  corresponding 
variations  of  rolling  moment  coefficient  versus  angle  of  attack  have  been  plotted 
in  Figures  24-26  for  JR  = 5.07,  8 and  12,  respectively.  The  figures  show  a hysteresis 
loop  for  rolling  moment  coefficient,  which  behaves  with  respect  to  lift  curve  slope 
and  aspect  ratio  in  a manner  very  similar  to  lift  coefficient. 


4.3.2  EFFECTS  OF  TAPER  RATIO,  SWEEP  ANGLE  AND  TWIST.  Effects  of 
taper  ratio,  sweepback  angle  and  twist  on  lift  hysteresis  were  obtained  for  the  At  =8 
case  with  section  lift  curves  #4,6  and  7.  Thus  figures  27-29  show  effects  of  varying 
t from  0.2  to  1.0  on  lift  hysteresis.  The  calculations  for  the  intermediate  value 
t = 0.495  are  the  same  as  those  presented  previously  in  Figure  22,  and  exhibit  a 
slightly  greater  tendency  for  lift  hysteresis  for  all  negative  lift  curve  slopes  than 
do  either  of  the  other  two  t cases.  This  may  be  true,  because  the  t = 0.495  taper 
ratio  provides  a more  uniform  spanwise  distribution  of  section  C^»  and  hence  causes 
the  initial  stall  to  be  delayed  to  higher  values  of  overall  lift  coefficient  when  the  pitch 
rate  is  positive. 


Lift  hysteresis  curves  with  a sweepback  angle  A = 20°  have  been  plotted  in  Figure  30 
for  section  lift  curves  #4,6  and  7.  Figure  30  may  be  compared  with  Figure  22  (a-c) 
to  evaluate  effects  of  the  increased  sweep  angle.  Although  these  figures  show  only 
a minor  sweep  effect,  it  should  be  noted  that  the  present  lifting  line  approach  may 
well  apply  only  for  low  sweep  angles,  because  of  the  assumption  in  the  present  analysis 
that  the  section  stalling  characteristics  remain  two-dimension.  The  introduction  of 
even  a moderate  degree  of  wing  sweep  is  known  to  modify  the  chordwise  pressure 
distribution  and  to  promote  spanwise  velocity  components,  thereby  compromising  the 
assumption  of  locally  two-dimensional  stall. 


Effects  of  wing  twist  angle  are  shown  in  Figures  31-33.  Each  figure  is  for  a different 
negative  lift  curve  slope  and  presents  results  for  the  cases  of  -3C  linear  wash  out 
(0  s ox  < -3°),  an  untwisted  wing  (^^-0),  and  a wing  with  positive  linear  twist 
(Osorj.  < +3°).  The  calculations  for  the  untwisted  case  are  the  same  as  presented 
previously  in  Figure  22(a-c).  As  is  apparent  from  Figures  31-33,  twist  has  only  a 
minor  influence  on  the  extent  of  life  hysteresis  for  lift  curves  #6  and  #7.  However, 
for  lift  curve  #4,  going  from  wash  out  ( a T = o -*  -3°)  to  wash  in  (a  T=  0 -*  +3°)  is  seen 
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to  increase  the  extent  of  stall  hysteresis.  This  may  be  due  to  the  more  uniform 
spanwise  distribution  of  section  CL  produced  by  wash  in  at  the  tip,  as  was  the  case 
with  the  effect  of  taper  ratio. 

It  is  recognized  that  the  effects  of  twist,  taper  ratio,  aspect  ratio,  and  sweep  back 
on  stall  hysteresis  are  interrelated.  The  particular  combination  of  wing  geometric 
parameters  which  provides  the  most  uniform  spanwise  distribution  of  section  C 
should  also  exhibit  a greater  degree  of  stall  hysteresis.  Because  of  the  correspondanoe 
between  the  occurrence  of  stall  hysteresis  and  zero  beta  rolling  and  yawing  moments 
(see  Section  4.3.3),  it  is  expected  that  similar  conclusions  may  be  made  with  regard 
to  the  effects  of  wing  planfom  parameters  on  the  zero-beta  moments. 

4.3.3  ZERO  BETA  ROLLING  (AND  YAWING)  MOMENTS.  A set  of  computations 
was  made  to  illustrate  the  occurence  of  zero-beta  rolling  and  yawing  moments  during 
the  pitch  sequence.  This  is  shown  in  Figures  34  and  35  for  the  T-2C  planform, 

(AS  = 5.07).  In  this  case  the  pitching  motion  was  initiated  at  Op  = 12°  with  the  wing 
symmetrical.  The  pitching  motion  was  halted  at  a p = 14. 4°  and  an  asymmetry 
A Or  = ±1°  was  applied  to  each  wing  plane  and  then  subsequently  removed,  in  order 
to  simulate  an  aileron  perturbation.  The  pitching  motion  was  then  resumed,  either 
positive  or  negative,  with  the  wing  again  symmetrical.  Figure  34  shows  the  lift 
variation  during  this  process  for  various  lift  curves.  The  corresponding  rolling 
moment  variations  are  shown  in  Figure  35.  The  theory  predicts  the  occurence  of 
zero-beta  rolling  and  yawing  moments  (not  plotted)  over  nearly  the  same  ap  zone 
as  that  for  which  lift  hysteresis  occurs.  Reducing  the  abruptness  of  the  stall 
(steepness  of  the  negative  lift  curve  slope)  is  seen  to  ameliorate  the  effect. 

A similar  set  of  computations  is  shown  in  Figures  36  and  37  for  an  A1  = 8 wing. 

Here  the  roll  asymmetry  was  introduced  and  removed  atop  = 12.4°,  since  the 
wing  was  already  completely  stalled  at  ap  = i4.4°.The  results  show  that  increasing 
aspect  ratio  from  5. 07  to  8. 0 considerably  reduces  the  zone  of  angle  of  attack  over 
which  zero- >eta  lateral  moments  occur,  especially  with  the  less  abrupt  section  lift 
characteristics. 


t 4.3.4  EFFECTS  OF  NUMBER  OF  ELEMENTS.  All  of  the  parametric  calculations 

* presented  so  far  were  made  with  14  vortex  elements  across  the  span,  in  order  to 

' conserve  computer  time.  It  is  recognized  that  requirements  for  improved 

! computational  accuracy  could  well  dictate  the  use  of  an  increased  number  N of  span- 

wise  elements  (the  existing  computer  program  is  limited  to  N = 20).  Figure  38  shows 
the  effect  of  varying  N from  10  to  20  on  lift  hysteresis,  for  the  section  characteristics 
with  the  most  abrupt  stall  (curve  #4)  and  with  the  most  gradual  stall  (curve  #8). 

I Although  changing  N appears  to  effect  the  instantaneous  Cl  value  somewhat,  especially 

? for  the  case  with  the  most  gradual  stall,  the  overall  stalling  characteristics  appear 

! virtually  uneffected  by  N.  It  is  recognized  that  further  studies  of  the  effect  of  N, 

; covering  a wider  range  well  beyond  N = 2o,  and  for  various  values  of  aspect  ratio 


and  taper  ratio,  are  needed  to  determine  an  optimum  value  for  the  computations. 

In  lieu  of  any  comprehensive  study  of  this  type,  the  value  N=14  used  in  the  majority 
of  the  computations  appears  reasonable. 

Most  of  the  parametric  computations  have  been  made  with  only  four  rows  of  vortex 
elements  in  the  wake  (M-4),  the  last  row  being  horseshoe  in  shape.  Calculations 
were  made  with  values  of  M as  high  as  40  (and  as  low  as  1).  The  value  M = 4 was 
set  as  a compromise  for  including  the  major  contribution  from  the  unsteady  wake 
(up  to  distances  of  3c  behind  the  airfoil),  while  reducing  computer  time  to  a minimum 

4.4  HYSTERESIS  CORRELATION 

An  attempt  was  made  to  correlate  the  results  of  the  various  computations  of  lift 
hysteresis  with  the  abruptness  of  the  wing  stall  and  with  planform  geometry.  One 
such  correlation  is  shown  in  Figure  39,  wherein  the  product  AC^  • A^h  has  been 
plotted  versus  57. 3b/7rAR,  Here  ACl^  is  the  maximum  C^  difference  in  the 
hysteresis  zone,  and  Aojj  is  the  width  of  the  zone.  The  product  ACl  • Affpj 
is  a measure  of  the  area  of  the  lift  hysteresis  loop.  The  value  b is  the  maximum 
negative  lift  curve  slope  (per  degree),  and  the  ratio  57.3  b/irAR  represents  the  ratio 
of  this  slope  to  the  slope  dCL/dag  for  a wing  with  constant  downwash  angle  crD 
across  the  span. 

It  is  readily  seen  from  Figure  1 that  multiple  solutions  (symmetric  case)  will  not 
occur  unless  b/(dCL/da s l,  corresponding  to  57.3  b/n  AR  £ 1 in  the  constant 
downwash  angle  case.  For  a wing  with  AR  - 4 and  a trilinear  lift  curve,  as  used 
in  the  examples  of  Section  2.4,  multiple  asymmetric  solutions  occur  when  b/'a 
(57.3  b)/(2~  AR/4)  £ 1.5,  (or  equivalently  57.3  b/v  AR  2 0.75),  as  seen  from  Figures 
7-9.  Thus,  the  parameter  57.3  bArAR  appears  reasonable  for  predicting  the 
occurence  of  multiple  asymmetric  as  well  as  of  symmetric  solutions. 

The  correlations  in  Figure  39  predict  the  occurrence  of  lift  hysteresis  wherever 
57. 3 b/ ir&>0. 18,  although  the  area  of  the  hysteresis  loop  is  seen  to  be  dependent 
on  other  planform  characteristics  as  well,  such  as  sweep  back  angle,  twist  and 
taper  ratio.  Thus,  the  method  predicts  the  occurrence  of  lift  hysteresis  for  wings 
with  negative  lift  curve  slopes  only  approximately  l/5  as  abrupt  as  thought  necessary 
based  on  the  idealized  constant  downwash  angle  criterion.  This  may  be  due  to  the 
ratio  dC^/do  being  much  less  than  it  £3/57.3  over  local  regions  of  the  wing  with 
stalled  or  partially  stalled  flow.  Additional  analysis  of  the  parametric  computations 
and  dynamic  wind  tunnel  testing  are  required  to  further  substantiate  the  correlation 
in  Figure  39. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


A nonlinear  lifting  line  theory,  with  an  unsteady  wake,  has  been  developed  and  has 
been  applied  to  the  prediction  of  forces  and  moments  on  aircraft  undergoing  a 
pitching  motion  through  the  stall  region.  The  procedure  is  an  outgrowth  of  the  steady 
state,  nonlinear  lifting  line  approaches  in  Reference  1-3,  which  predicted  the  occur- 
rence of  multiple  solutions  to  the  lifting  line  equations  under  certain  conditions.  The 
present  theory  predicts  the  occurrence  of  lift  hysteresis  and  of  rolling  and  yawing 
moments,  even  at  zero  yaw  angle,  for  wings  where  section  characteristics  exhibit 
a negative  lift  curve  slope  of  sufficient  steepness  subsequent  to  stall.  It  is  anticipated 
that  the  theory  (and  parametric  computations  presented  herein)  may  serve  as  a guide 
in  predicting  the  susceptibility  of  various  aircraft  designs  of  moderate  to  high  aspect 
ratio  and  low  sweep  angle  to  adverse  stall  and  departure  characteristics. 

5.1  CONCLUSIONS 

Conclusions  relating  to  the  numerical  computations  and  to  the  assumptions  in  the 
theory  are  summarized  below. 

(i)  The  method  has  been  programmed  for  solution  on  the  CDC  Cyber  70 
computer.  Computer  solutions  are  simple  to  execute  and  reasonably 
economical,  requiring  about  one-minute  of  micro-processor  time  for  a 
complete  pitch  sweep  and  hysteresis  loop. 

( ii)  Computations  by  the  present  method  agree  reasonably  well  with  previous 
more  exact,  unsteady,  linear,  two-dimensional  and  three-dimensional 
theories  down  to  aspect  ratios  of  3 and  less  for  wings  of  zero  sweep. 

(iii)  The  present  method  compares  well  with  steady  state  T-2C  test  data  (aspect 
ratio  5.  07),  both  for  lift  and  rolling  moment.  No  lift  hysteresis  effect 
appeared  in  the  test  data,  but  both  test  data  and  theory  showed  the  occurrence 
of  zero-beta  rolling  and  yawing  moments  in  the  stall  region. 

(iv)  The  present  theory  predicts  the  occurence  of  lift  hysteresis  (and  of 

zero  beta  rolling  and  yawing  moments)  when  the  parameter  57. 3b/71, At  > 0. 18 
(-bis  the  maximum  negative  lift  curve  slope  of  the  two-dimensional  airfoil 
sections  comprising  the  wing).  The  angle  of  attack  range  over  which  these 
effects  occur,  and  their  severity  is  found  to  Increase  with  increasing  b 
and  with  decreasing  aspect  ratio.  Other  planform  parameters,  such  as  sweep- 
back  angle, twist  and  taper  ratio  have  also  been  shown  to  affect  the  extent 
of  these  adverse  stalling  characteristics  to  a lesser  degree. 
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(v)  The  converged  solutions  obtained  at  angles  of  attack  for  which  multiple 
roots  exist  are  non  unique  and  depend  upon  the  starting  solutions  used  in 
the  integration  procedure. 

(vi)  The  iteration  procedure  used  in  the  present  method  employed  a special 
logic,  which  used  either  the  converged  solution  from  the  previous  time 
step  as  the  initial  iterate  (unstalled)  case,  or  set  the  initial  iterate  for 
downwash  angle  equal  to  zero  over  a wing  panel  (stalled  case).  This 
special  logic  was  somewhat  empirical,  and  led  to  the  elimination  of  solu- 
tions with  saw-tooth  stalling  patterns.  However,  the  special  logic 
delayed  stall  during  pitch-up  motions,  and  postponed  destall  during  pitch- 
down  motions,  thereby  maximizing  the  zone  over  which  lift  hysteresis 
(and  zero  beta  rolling  and  yawing  moments)  occur. 

(vli)  The  formulation  does  not  Include  geometric  effects  of  sweep  angle  and 
dihedral  angle  on  the  effective  section  angle  of  attack.  These  effects 
should  be  included  to  produce  satisfactory  correlation  with  test  data, 
particularly  for  the  s idesllp  case. 


5 2 RECOMMENDATIONS 

Recommendations  for  improving  the  present  method,  by  removing  several  of  the 
restrictive  assumptions  noted  above,  by  conducting  additional  comparisons  with  test 
data,  and  by  generalizing  the  procedure  for  a wider  class  of  aircraft  configurations 
(e.g.  lower  aspect  ratio  fighters)  are  given  below. 

(i)  A less  empirical  approach  is  required  for  selecting  a particular  loading 
when  multiple  solutions  exist.  This  would  lead  to  an  improved  prediction 

of  the  angle  of  attack  zone  over  which  lift  hysteresis  and  attendant  zero-beta 
lateral  moments  may  exist.  One  approach  may  be  the  introduction  of  rate 
dependent  equations  for  the  vortlcity  on  the  wing.  Such  a formulation  may 
require  chordwise  relocation  of  the  control  point  and/or  modification  of  the 
time  stepping  procedure  (see  also  iv  below). 

(ii)  A more  definative  evaluation  of  the  method  should  be  carried  out  by  conducting 
wind  tunnel  tests  with  a model  which  may  be  pitched  and  yawed  continuously 
while  obtaining  instantaneous  readings  of  forces  and  moments.  It  is 
recommended  that  the  model  consist  of  a circular  fuselage  with  replaceable 
straight  wings  of  various  aspect  ratio  and  of  known  two-dimensional  section 
characteristics,  in  order  to  determine  limits  of  applicability  of  the  method. 
Pressure  data  should  be  obtained  over  the  wing,  so  that  transient  span  and 
chordwise  load  distributions  may  be  obtained  and  compared  with  those 
predicted  by  the  theory. 
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(Ill)  The  formulation  should  be  extended  to  treat  yaw  and  dihedral  effects  directly, 
rather  than  by  simulating  a yawed  wing  by  a skewed  wing.  In  addition,  a 
capability  for  time  dependent  yawing  motions  should  be  incorporated  into  the 
program  In  order  to  provid  2 for  coupling  the  lifting  line  formulation  with  the 
equations  of  motion  for  an  aircraft. 

(iv)  In  order  to  extend  the  method  to  lower  aspect  ratio  fighter  type  wings,  it  Is 
recommended  that  the  lifting  line  formulation  be  expanded  into  an  unsteady, 
nonlinear,  lifting  surface  theory.  The  new  formulation  would  employ  at 
least  two  spanwise  lifting  lines,  and  control  points  for  satisfying  flow  tangency 
would  be  placed  at  two  or  more  chordwise  positions.  The  expanded  theory 
would  also  include  unsteady  effects  on  chordwise  load  distribution,  and 
might  provide  insight  into  the  mechanism  of  dynamic  stall.  The  later  effects 
are  excluded  from  the  current  lifting  line  formulation,  which  assumes  that 
chordwise  load  distributions  are  always  steadystate. 
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SECTION  LIFT  COEF,  C,  SECTION  LIFT  COE F. 


ANGLE  OF  ATTACK,  Cl 
A.  TYPICAL  SECTION  LIFT  CURVE 


ANGLE  OF  ATTACK, a 

B.  TRI-LINEAR  SECTION  LIFT  CURVE 
USED  FOR  SIMPLIFIED  CALCULATIONS 


Figure 

<T< 

y 


Nonlinear  Two-Dimensional  Section  Lift  Curves  Which  Yield 
Multiple  Solutions  of  the  Lifting  Line  Equations 
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Figure  3.  Vortex  Elements  for  Unstesoy  Wake  Formulation 
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Figure  5.  Two-Dimensional  Representations  of  Unsteady  Wake 
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AC,  (t)  / 27: Aa 
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0.4 


® PRESENT  METHOD 
HAND  COMPUTATIONS 
n = 1 (STEP  SIZE  EQUAL 
TO  ONE  CHORD  LENGTH) 

f = 0.25: 


0.2 


DISTANCE  TRAVELED  BY  WAKE  BEHIND 
TRAILING  EDGE  IN  CHORD  LENGTHS 
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Figure  6.  Comparison  of  Discrete  Vortex  Wake  Model  with  Wagner  Function 
(Continuous  Vortex  Wake)  for  Two-Dimensional  Flow 


Figure  9.  Asymmetric  Loadings  for  Wings  with  One  Panel  Unstalled  (dc^/da 
and  the  Other  Fully  Stalled  (dC^/da  =0) 


MAX 


TABULATION  OF  SYMMETRICAL  ANO  ASYMMETRICAL  LOADINGS  FOR 
2 PANEL  SOLUTION,  AR=4  WING  WITH  b/a=4,  AC. /C,  =0.7,  ANDO- /a,=1.3. 
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Figure  10.  Example  of  the  Occurance  of  Six  Possible  Loadings  with  a Trl-Llnear 
Lift  Curve 
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Figure  11.  Predicted  Stability  and  Convergence  Limits  for  Iteration  Procedure 
with  a Single  Vortex  Element 
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RECTANGULAR  WING 
ASPECT  RATIO  6 

SECTION  LIFT  CURVE  SLOPE  a = 2r? 


STEADY 

STATE 


J I I I L 

1 2 3 4 5 

DISTANCE  TRAVELED,  Ax/c 


Figure  12.  Effect  of  a Step  Change  in  Angle  of  Attack  on  Lift,  AR  - 6. 
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SOLID  LINES  BY  CURRENT  METHOD 
RECTANGULAR  WING 
N=  14,  M=  20,  n=  1 
SECTION  LIFT  CURVE  SLOPE  a=  2rr 


0«- 

0 


Figure 
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14. 


2 4 6 8 10 

DISTANCE  TRAVELED,  Ax/c 

Effect  of  Aspect  Ratio  on  Lift  Following  a Step  Change  in  Angle  of 
Attack 
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Figure  16.  Planform  Geometry  and  Paneling  Used  for  T-2C  Computations  (N  = 14) 


MOMENT  COEFFICIENT,  CM  UFT  COEFFICIENT,  C 


CASE  QlNTENOEO  TO  REPRESENT  NACA  64i  A212  AIRFOIL  Re  = 3 X 106 
CASE  0INTENOED  TO  REPRESENT  NACA  641  A212  AIRFOIL  Re=  6 X 
CASE  0 INTEND  ED  TO  REPRESENT  NACA  64i  A212  AIRFOIL  Re>6X  10^ 


RATA  SETS  11, 12,  AND  13  NOTED  BY  As  FOR  PRE  STALL 
0's  & SAME  AS  1,  2 & 3 FOR  CORRESPONDING  POST  STALL 
(U  THESE  REPRESENT  DATA  FOR -3°  T.E.  DEFLECTION. 


ANGLE  OF  ATTACK,  ap  'degrees! 


A.  LIFT  COEFFICIENT 


for  CASE 


ANGLE  OF  ATTACK, ap  (degrees) 

B.  PITCHING  MOMENT  COEFFICIENT  ABOUT  25%  CHORD  POINT 


Figure  17.  Section  Characteristics  of  NACA  64  A212  and  T-2C  Airfoils 


ANGIE 

C.  P30 

igure  17.  Section  Characteristics  of I 


MOMENT  COEFFICIENT,  LIFT  COEFFICIENT,  C 


CASE  4 REPRESENTS  Re  = 3 X 106,  NACA  64-209 
CASE  5 REPRESENTS  Re  = 6 X 106,  NACA  64-209 
CASES  G-8  REPRESENT  LESS  ABRUPT  STALLS 


S 12  16  20 

2 

ANGLE  OF  ATTACK,  ap  (degrees) 
A LIFT  COEFFICIENT 

ANGLE  OF  ATTACK, ap  (degrees) 

B.  PITCHING  MOMENT  COEFFICIENT  A80U  T 25%  CHORD  POINT 


Figure  18.  Nonlinear  Airfoil  Section  Characteristics  Used  for  Parametric 
Calculations 
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0 I.  I l — ,1,  ..  ■' 

0 4 8 12  16  20  24 

ANGLE  OF  ATTACK,  ap  (degrees) 


C.  PROFILE  DRAG  COEFFICIENT 

Figure  18.  No  .linear  Airfci!  Section  Characteristics  Usee!  for  Parametric 
Calculations  (Concluded) 
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1.5 


O TEST  DATA 


PRESENT  METHOD  WITH  SECTION 


XI 1 


A.  ZERO  YAW  ANGLE 


T-2C  DATA  FROM  REF.  23. 
N.c  = 0.20,  Re-  = 4 X 106 
GRiTOFF 


B.  ±10°  YAW  ANGLE 

Figure  19.  Comparison  of  Present  Method  with  T-2C  Wind  Tunnel  Data  for  Lift 
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, ROLLING  MOMENT  COEFFICIENT  ±ACVR  (degrees) 


TEST  DATA  FROM  REF.  23 
Re  =*4x  i0*>  BASED  ON  c 
O GRIT  OFF,  AGRITON 

CALCULATED  BY  PRESENT  METHOD 
SHOWN  BY  LINE  WITH  ARROWS 
N = 20,  M=  1,  SECTION  DATA 
FROM  CURVE  NO.  11  OF  FIGURE  17 


Figure  20.  Compailsonof  Present  Method  with  T-2C  Wind  Tunnel  Data  for 
Rolling  Moment  at  Zero  Yaw  Angle 


ap  (degress) 

a NEGATIVE  LIFT  CURVE  SLOPE  « ai25PER  DEGREE 

Figure  21.  Effect  of  Negative  Lift  Curve  Slope  on  Lift  Hysteresis,  Aspect 
Ratio  5.07 


2D  LIFT  CURVE  NO.  4 
FROM  FIGURE  18^ 


A.  NEGATIVE  LIFT  CURVE  SLOPE  = (L97  PER  DEGREE 


AR  = 8.0,  t = 0.495,  Z27” 

0°,  aT=  0°,  7=  0°,  RB=  0.119  h/2 
PITCH  RATE  8 DEG/SEC 
SIDESLIP  = 0,  AaR=  0 
N = 14,  M = 4 


2D  LIFT  CURVE  NO.  6 
FROM  FIGURE  ^ 


ap  (degree) 

1 NEGATIVE  LIFT  CURVE  SLOPE  = 0.125 PER  DEGREE 

Figure  22.  Effect  of  Negative  Lift  Curve  Slope  on  Lift  Hysteresis,  Aspect 
Ratio  8. 0 
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20  LIFT  CURVE  NO.  7 
FROM  FIGURE  18^ 


C.  NEGATIVE  LIFT  CURVE  SLOPE  = 0.075 


NOTATIONS  ON  FIGURES  A & B 


AR  = 007,  t-  01495,  A=  2.27 
if  0°,  ay*  0°,  y=  0°,  Rb  = a 119  b/2 
PITCH  RATE  8 DEG/SEC 
SIDESLIP- 10  DEG,  AaR=0° 

N = 14,  M « 4 


2D  LIFT  CURVE  NO.  4 
pphu  nniiRr  in 

rnUM  riuUnL  lo 

r*V4 TA  A A , 


A.  NEGATIVE  LIFT  CURVE  SLOPE  = 097  PER  DEGREE 


2D  LIFT  CURVE  NO.  6 
FROM  FIGURE  18 


0 NEGATIVE  LIFT  CURVE  SLOPE  = 0125  PER  DEGREE 


20  LIFT  CURVE  NO.  7 

- a.  . r 

FROM  FIGURE  18 

C.  NEGATIVE  LIFT  CURVE  SLOPE  = 0075  PER  DEGREE 


2D  LIFT  CURVE  NO.  8 
FROM  FIGURE  18 


D.  NEGATIVE  LIFT  CURVE  - 0040  PER  DEGREE 

Figure  24.  Effect  of  Negative  Lift  Curve  Slope  on  Rolling  Moment  Hysteresis  at  10 
Sideslip  Angle, Aspect  Ratio  5.07 


AR  = 12.0,  r = 0.495,  A=  2.27° 

^ 0°,  aT=  0°,  7=  0°,  Rb=  0.119  b/2 
PITCH  RATE  8 DEG/SEC 
SIDE  SLIP-'  10  DEG 
N - 14,  M ” 4 


B.  NEGATIVE  LIFT  CURVE  SLOPE  = ai25  PER  DEGREE 


D.  NEGATIVE  LIFT  CURVE  SLOPE  = G040  PER  DEGREE 


Figure  2G.  Effect  of  Negative  Lift  Curve  Slope  on  Rolling  Moment  Hysteresis  at  10° 
Sideslip  Angle,  Aspect  Ratio  12.0 
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AR  = &0,  A=  227°,  i = 0 

aT=  0°,  7=  0°,  Rb=  0.119  b/2 

PITCH  RATE  = 8 DEG/SEC,  b=  a 125  PER  DEG 


Figure  28.  Effect  of  Taper  Ratio  on  Lift  Hysteresis,  Aspect  Ratio  8. 0,  Lift  Curve 
#6  from  Figure  18 
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ap  (dtgrees) 
a TAPER  RATIO  1.00 

Figure  29.  Effect  of  Taper  Ratio  on  Lift  Hysteresis,  Aspect  Ratio  8.0,  Lift  Curve 
#7  from  Figure  18 


AR  = 5.07, r = 0.495,  A=  2.27  SIDESLIP 

i = 0°,  aj  = 0°,  7=  0°,  Rg=  0.119  b/2  Aa„  = 0° 
PITCH  RATE  = ±8  OEG/SEC  N=14.M  = 


SOLID  SYMBOLS  WITH 
Aa0  - ±1°  OH  L/R  WING 


2D  LIFT  CURVE  NC.  4 
FROM  FIGURE  IS 


A.  NEGATIVE  LIFT  CURVE  SLOPE  = (L97  PER  DEGREE 


v v v v v 

X V 

20  LIFT  CURVE  NO.  6 

FRGfA  FIGURE  18 

* 

~ 1 ! 

I i 

8.  NEGATIVE  LIFT  CURVE  SLOPE  = 0.125  PER  DEGREE 


2D  LIFT  CURVE  NO.  7 
FROM  FIGURE  18 


C.  NEGATIVE  LIFT  CURVE  SLOPE  = 0075  PER  DEGREE 


20  LIFT  CURVE  NO.  8 
FROM  FIGURE  18 


0^  (degrees) 

D.  NEGATIVE  LIFT  CURVE  SLOPE  = 0.040  PER  DEGREE 


Figure  35.  Effect  of  Roll  Perrurbauoc  at  c*p  14.4  Degrees  on  Zero  Beta  Rolling 
Moment,  Aspect  Ratio  5.07 


AR  = 8.0,  r = 0.495,  A=  227 
i = 0o,aT=0°,r=0°,  RB=  ai19  b/2 
PITCH  RATE  =±8  DEG/SEC 
SIDESLIP  = 0° 

Aafi  = 0°  (OPEN  SYMBOLS) 

N = 14,  M = 4 


SOLID  SYMBOLS  WITH 
AaR  = ±1°  ON  L/R  WING 
PANELS  AT  a = 124° 


20  LIFT  CURVE  NO.  4 
FROM  FIGURE  18 


A.  NEGATIVE  LIFT  CURVE  SLOPE  = (197  PER  DEGREE 

I I I 


2D  LIFT  CURVE  NO.  6 
FROM  FIGURE  38 


B.  NEGATIVE  LIFT  CURVE  SLOPE  = M25  PER  DEGREE 


)O0O0° 


20  LIFT  CURVE  NO.  7 
FROM  FIGURE  18 


C.  NEGATIVE  LIFT  CURVE  SLOPE  = Q.075  PER  DEGREE 

till 


20  LIFT  CURVE  NO.  8 
FROM  FIGURE  18 


0.  NEGATIVE  LIFT  CURVE  SLOPE  = (L04  PER  DEGREE 
I I I 


ap  (degree) 


Figure  36.  Effect  of  Roll  Perturbation  atQ'p  = 12.4  Degrees  on  Lift  Hysteresis 
Aspect  Ratio  8. 0 ai 


AR  = 8.0,  r = 0.495,  A=  2.27° 
i=  0°,  aT=  0°,  y=  0°,  Rg  = Oil  19  b/2 
PITCH  RATE  = ±8  DEG/SEC 
AaR  = 0°  (OPEN  SYMBOLS) 

N = 14,  M = 4 


SOLID  SYMBOLS  WITH 
AaR  = ±1°  ON  L/R  WING 
PANELS  AT  dp  = 12,4° 
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6 8 10  12  14  16 

ap  (degrees) 


0.  NEGATIVE  LIFT  CURVE  SLOPE  = 0.040  PER  DEGREE 


Figure  37.  Effect  of  Roll  Perturbation  at  Op  = 12.4  Degrees  on  Zero  Beta  Rolling 
Moment,  Aspect  Ratio  8.0 
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2P  LIFT  CURVE  NO.  8 
FROM  FIGURE  18 


N = 20 
N=  12 

N=  10  OMITTED 
FOR  CLARITY 


